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DISPERSIVE AND DIFFUSIVE-DISPERSIVE SHOCK WAVES FOR 
NON-CONVEX CONSERVATION LAWS* * * § 
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Abstract. We consider two physically and mathematically distinct regularization mechanisms 
of scalar hyperbolic conservation laws. When the flux is convex, the combination of diffusion and 
dispersion are known to give rise to monotonic and oscillatory traveling waves that approximate shock 
waves. The zero-diffusion limits of these traveling waves are dynamically expanding dispersive shock 
waves (DSWs). A richer set of wave solutions can be found when the flux is non-convex. This review 
compares the structure of solutions of Riemann problems for a conservation law with non-convex, 
cubic flux regularized by two different mechanisms: 1) dispersion in the modified Korteweg-de Vries 
(mKdV) equation; and 2) a combination of diffusion and dispersion in the mKdV-Burgers equation. 
In the first case, the possible dynamics involve two qualitatively different types of DSWs, rarefaction 
waves (RWs) and kinks (monotonic fronts). In the second case, in addition to RWs, there are 
traveling wave solutions approximating both classical (Lax) and non-classical (undercompressive) 
shock waves. Despite the singular nature of the zero-diffusion limit and rather differing analytical 
approaches employed in the descriptions of dispersive and diffusive-dispersive regularization, the 
resulting comparison of the two cases reveals a number of striking parallels. In contrast to the 
case of convex flux, the mKdVB to mKdV mapping is not one-to-one. The mKdV kink solution 
is identified as an undercompressive DSW. Other prominent features, such as shock-rarefactions, 
also find their purely dispersive counterparts involving special contact DSWs, which exhibit features 
analogous to contact discontinuities. This review describes an important link between two major 
areas of applied mathematics, hyperbolic conservation laws and nonlinear dispersive waves. 


1. Introduction. Shock waves represent one of the most recognizable features 
of nonlinear wave systems. In their simplest form, shock waves can be modeled as 
discontinuous, weak solutions of scalar conservation laws in one space dimension x 
and time t 


u t + f(u) x = 0. (1.1) 

When the flux function / is convex, shock admissibility rests upon an entropy condi¬ 
tion that uniquely selects the physically relevant weak solution [64]. The existence of 
entropy presupposes some irreversible dissipative mechanism, e.g., diffusion, in the un¬ 
derlying physical medium. A simple example of a diffusively regularized conservation 
law is Burgers’ equation 


Ut 2 — VUxxi ( 1 - 2 ) 

where v > 0 is a constant representing the diffusion strength. 

However, in conservative physical systems whose dominant shock regularizing 
mechanism is dispersion, there is no corresponding notion of entropy. Discontinuous, 
weak solutions of (1.1) no longer accurately model wavebreaking dynamics. An alter¬ 
native theory of dispersive shock waves (DSWs) [45] seeks to examine the asymptotic 
structure of solutions of dispersive regularizations of (1.1), the simplest example being 
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the Korteweg-de Vries (KdV) equation 

Ut V 2 = f-^ U-XXX , (1*3) 

in which /r 7 ^ 0 is a constant. 

The Burgers’ (1.2) and KdV (1.3) equations represent two distinctly different 
ways to regularize the conservation law (1.1) with the convex flux f(u) = \u 2 . When 
the flux in ( 1 . 1 ) is non-convex, new wave features emerge including composite wave 
solutions and undercompressive shocks that do not satisfy the Lax entropy condition 
[ 68 ]. Many authors have approached the issue of shock wave admissibility when the 
flux is non-convex by introducing criteria based on mechanics or physics, notably 
through kinetic relations and nucleation conditions [1, 69]. An additional approach 
includes regularization of the equation by higher order terms that are dissipative and 
dispersive, representing more of the physics than is contained in the conservation law 
(1.1) alone [ 8 , 68 ]. The analysis of shock dynamics in dispersive equations is in many 
respects quite different from the theory of conservation laws, as exemplified by the 
famous book of Whitham [97], the pioneering work of Gurevich and Pitaevskii [40] and 
the series of papers of Lax and Levermore [65]. Much less is known about dispersive 
shock waves in equations with non-convex flux. 

The aim of this paper is to review and compare results of these two regularization 
approaches for shock wave solutions of (1.1). The comparison is well known for the 
case of convex flux functions (see, e.g., [46]), but it reveals a number of new nontrivial 
parallels and contrasts if / is non-convex. This comparison serves two purposes: 1) 
it identifies the dispersive theory as a singular limit of the diffusive-dispersive theory, 
2 ) it serves as a bridge between dispersive nonlinear waves and hyperbolic conser¬ 
vation laws, two major areas of applied mathematics. As far as possible, dispersive 
shock waves are described in the language of hyperbolic conservation laws (charac¬ 
teristics, admissibility criteria, non-classical, undercompressive shocks) and, similarly, 
diffusive-dispersive shocks are explored in the language of dispersive shock theory 
(shock polarity and orientation). Most of the results presented here exist in the lit¬ 
erature. The purpose of this review is to distill them into a cohesive narrative that 
contrasts and unifies two areas of applied mathematics. 

We focus on the modified KdV-Burgers (rnKdVB) equation 

T (rt )a; = VUxx + l^'^XXX') (1'4) 

in which v > 0 and /i / 0 are parameters. When v = 0, the equation is purely 
dispersive, and we refer to the equation using the abbreviation rnKdV. For v > 0, the 
equation has both diffusion and dispersion and we refer to it as the rnKdVB equation. 
The modified equation is considered because the combination of non-convex flux with 
dispersion yields a rich collection of wave solutions. Furthermore, eq. (1.4) can be 
viewed as a universal model of diffusive-dispersive nonlinear wave dynamics under a 
certain criticality condition. We demonstrate this by deriving (1.4) from a class of 
diffusive-dispersive Euler equations. 

The applied significance of purely dispersive conservation law regularization is 
twofold. Firstly, DSW dynamics describe the initial stage of shock development in 
media with both regularization mechanisms present (subject to an appropriate scale 
separation of diffusive and dispersive effects). However, the most interesting and 
rich applications of DSW theory arise in media whose regularization dynamics are 
dominated by dispersion alone. Prominent examples of such media are Bose-Einstein 
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condensates and nonlinear optical media (see, e.g., [46], [84], [95] for relevant ex¬ 
perimental and theoretical DSW studies), although conservative nonlinear-dispersive 
dynamics are also common in water waves and the atmosphere, where DSW-like 
structures, termed undular bores, have been observed [90], [13]. 

1.1. Scalar conservation laws. The study of conservation laws goes back to 
Riemann’s analysis of the classic shock tube problem in which a membrane separating 
gas under different pressures is ruptured, generating a shock wave and a rarefaction 
wave propagating in opposite directions away from the membrane. This problem 
is simply represented by the nonlinear system (2.1) of isentropic gas dynamics (see 
§2), in the diffusion and dispersion-free limit (y = 0, D = 0). The study of such 
hyperbolic systems in one space dimension and time has progressed rapidly since the 
1950’s, with early work of Lax, Oleinik and others establishing existence of solutions 
of Riemann problems, shock formation and long-time behavior of solutions, for initial 
data with small variation. The seminal paper of Glimm [31] injected profound new 
ideas that were used to prove existence of solutions of initial value problems. As in 
Glimm’s paper, the method now referred to as wave front tracking [10, 15] relies on 
solutions of Riemann problems to approximate weak solutions of hyperbolic systems. 
Development of the theory in more than one space dimension has been more difficult, 
but papers of Majda [73, 74] have led the way to further refinements [37]. 

The theory of scalar conservation laws is much more complete than for systems. 
For convex fluxes, the theory of initial value problems in one space dimension is 
classical [63]. The existence and uniqueness results of Kruzhkov [60] in any number 
of dimensions rely on entropy conditions. 

A function u = u(x,t) in L°°(R. x R + ) is a weak solution of the initial value 
problem 


ut + f(u) x = 0, t > 0, x G R. u(x, 0) = g{x), i£l, 
if for every test function </> G C 1 (K x R + ) with compact support, we have 

nOO nO O 

I (ucf>t + f(u)(j) x ) dxdt + / g(x)<f>(x 1 0) dx = 0. (1-5) 

—oo J — oo 

Here, R + = [0, oo), so that the support of <p can include a portion of the initial line 
E x [i = 0], and <p can be non-zero there. When / is smooth and uniformly convex 
(i.e., there is 9 > 0 such that f"(u ) > 6 for all u € R), the initial value problem has 
a unique weak solution whose shocks satisfy the Rankine-Hugoniot conditions (see 
eq. (3.1) below), but also an entropy condition of Oleinik [79, 80]: there is a constant 
C > 0 with the property that 



u(x + z, t) — u(x, t) < (7(1 + l/t)z for almost all i£l, z > 0, t > 0. 

A readable account of the proof of this result is given in the text of Evans [27]. 
DiPerna [17] proved asymptotic results for large t for entropy solutions of systems of 
equations, following earlier analysis of Glimm and Lax [32]. The DiPerna results for 
scalar equations are proved succinctly in Evans’ text. 

Another analysis approach is to modify the equations with terms that make initial 
value problems well posed in suitable function spaces. Typically this regularization 
will be controlled by a small parameter, e, say; when e = 0, the equations revert to the 
original hyperbolic equations. This is generally a singular limit and the convergence 
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is to a function that is less smooth than the approximating solutions. The analysis 
therefore proceeds in two steps, the first being to show convergence in a suitable 
topology, the second to show that the limiting function is indeed a weak solution. 
Additional steps may be needed to show uniqueness of the limit, and to show that 
the solution satisfies appropriate side conditions such as entropy conditions. Even for 
scalar equations, the nature of the regularization, and convexity or lack of convexity 
of the flux, play crucial roles in this approach. 

In this paper, we focus on two kinds of regularization, involving either purely 
dispersive terms, or a combination of dispersive and diffusive terms. For scalar equa¬ 
tions with convex flux, these two cases yield strikingly different approximations of 
shock waves; the differences become even more profound as the flux becomes non- 
convex. However, there is also considerable correspondence between the structure of 
regularized shocks in the purely dispersive case and the structure in the case of mixed 
dispersive-diffusive regularization. 

1.2. DSWs and modulation equations. A DSW manifests as a nonlinear, 
expanding wavetrain connecting two disparate states of slowly varying (or constant) 
flow. When viewed locally, i.e., over a small region of space and time, DSWs dis¬ 
play periodic, or quasiperiodic structure, forming due to the interplay between non¬ 
linear and dispersive effects. However, over a larger region covering multiple wave 
oscillations, the DSW wavetrain reveals slow modulation of the wave’s parameters 
(amplitude, frequency, mean etc.), and this modulation itself behaves as a nonlinear 
hyperbolic wave. This kind of “dispersive-hyperbolic” duality of modulated waves is 
familiar from linear wave theory, but DSWs prominently display it over a full range 
of nonlinearity, from a weakly nonlinear regime to solitary waves realized as an in¬ 
tegral part of the modulated wave train (see figure 1.1). The striking manifestation 
of nonlinearity in a DSW is that it is characterized by at least two distinct speeds of 
propagation, those of its leading and trailing edges. In contrast, energy transport by 
a linear wave packet is described by a single group velocity. The unsteady dynamics 
of DSWs have far-reaching physical and mathematical implications, among which are 
the principal inapplicability of the classical Rankine-Hugoniot relations and insepara¬ 
bility of the macroscopic DSW dynamics from the analysis of its nonlinear oscillatory 
structure. 

The mathematical description of DSWs involves a synthesis of methods from hy¬ 
perbolic quasi-linear systems, asymptotics, and soliton theory. One of the principal 
tools is nonlinear wave modulation theory, often referred to as Whitham averaging 
(Whitham 1965). Whitham theory was first applied to DSWs in the framework of the 
KdV equation by Gurevich and Pitaevskii [40]. The Gurevich-Pitaevskii approach 
rests upon the fundamental assumption that the DSW can be asymptotically repre¬ 
sented as a slowly modulated periodic traveling wave solution of the original nonlinear 
dispersive equation, e.g., the KdV equation (1.3), where the spatio-temporal modula¬ 
tion scale (variations of the wave amplitude, wavelength, mean etc.) is assumed to be 
much greater than the dispersion length, the wavelength scale of the carrier, locally 
periodic traveling wave (see figure 1.1). This scale separation enables one to effectively 
split the DSW description problem into two separate tasks of different complexity: 
the relatively easy problem of the periodic traveling wave solution description and the 
harder problem of finding an appropriate modulation which provides a match between 
the modulated wave train (the DSW) with the smooth, slowly varying external flow. 

Following Whitham [97], the equations describing slow modulations of periodic 
nonlinear waves can be obtained by averaging conservation laws of the original dis- 
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small amplitude waves 


leading edge 
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Fig. 1.1. The structure of a DSW with negative dispersion. The amplitude and wave number 
within the modulated wave are denoted a(x, t), k(x, t), respectively, with limits indicated at the leading 
and trailing parts of the DSW. The wave speeds s+, s~ refer respectively to the leading solitary wave, 
and to the trailing waves. The speed S- may be positive or negative. 


persive equations over the family of periodic traveling wave solutions. Given an N-th 
order nonlinear evolution equation, implementation of the Whitham method requires 
the existence of a iV-parameter family of periodic traveling wave solutions u) 
with phase 9 = x — Ut , phase velocity U( u), and parameters u £ K^. Additionally, 
the evolution equation must admit at least N—l conserved densities Vi[<f>] and fluxes 
Qi[(^>], i = 1,... N — 1 corresponding to the local conservation laws 

tV + 1 a = °- i = 1 ’ - w - L < L6 > 

We note that in the majority of known, physically relevant dispersive-hydrodynamic 
systems, the number of parameters characterizing periodic waves is N = 3 or N = 4, 
although systems yielding N > 4 can arise in some applications (see [24]). For the KdV 
and mKdV equations, N = 3 while for the more general case of Eulerian dispersive 
hydrodynamics considered in Section 2, N = 4. 

Assuming slow evolution of the wave’s parameters u = u(x, t) on spatio-temporal 
scales much larger than the wave’s wavelength L( u) and period T(u) = L/U , the 
conservation laws are then averaged over a wavelength resulting in the modulation 
equations 

(zjf 7?i ^ (0;u)]d0 ) + [ Q ^ u ')\ d ^j = 0 ’ * = 1, ■ • •, AT — 1. (1.7) 

The N—l averaged local conservation equations (1.7) are completed by the 
addition of the conservation of waves 

kt + = 0 , (1-8) 

where k(u) = 2t t/L and w(u) = kU = 2n/T are the wavenumber and frequency of 
the nonlinear traveling wave, respectively. Equation (1.8) represents a consistency 
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condition for the application of modulation theory. Note that the system (1.7), (1.8), 
unlike (1.6), is a system of non-dispersive conservation laws, which, assuming non¬ 
vanishing of the relevant Jacobians, can be represented in the standard form of a 
system of first order quasilinear equations 


u t + A(u)u x = 0. (1.9) 

The matrix A(u) encodes information about both the nonlinear and dispersion prop¬ 
erties of the original evolution equation. An alternative version of the Whitham 
method, leading to the same set of modulation equations (1.9), employs an averaged 
variational principle [98]. 

In some special but important cases, the system (1.9) can be reduced to diagonal 
form 

9 ^t +Vl{r) % z = °’ * = ( L1 °) 

where r.; (u-[ ,..., Ujv) are Riemann invariants, and the characteristic speeds V)(r) are 
eigenvalues of the Jacobian matrix A(u(r)). The existence of Riemann invariants 
dramatically simplifies the modulation equation analysis but is only guaranteed if 
N < 2. Remarkably, Whitham found Riemann invariants for the modulation system 
associated with periodic solutions of the KdV equation (1.3) where N = 3 [97]. A 
similar set of Riemann invariants was found in [18] for the rnKdV equation connected 
with the KdV equation by the Miura transform [76]. Later it was shown that the 
availability of Riemann invariants for the KdV-Whitham system is intimately linked 
to the integrable structure of the KdV equation via the inverse spectral transform 
[28], [58]. Subsequently, Riemann invariants were found for other modulation systems 
associated with integrable equations such as the nonlinear Sclrrodinger equation [29], 
[81], the Kaup-Boussinesq system [23] and others. We note that the notion of hyper- 
bolicity plays an important role in the analysis of modulation systems: for strictly 
hyperbolic systems, the mapping u i—> r is invertible and the characteristic speeds are 
distinct V\ < V 2 < ■ ■ ■ < Vn- We also note a useful representation of the characteristic 
speeds [39], [62], [51] 


Vi = 


1 - 


A A 
H dri 


U{ r), 


( 1 . 11 ) 


which follows from the consideration of the wave conservation equation (1.8) as a 
consequence of the Riemann system (1.10). 

The Whitham method can be viewed as a field-theoretic analogue of the Krylov- 
Bogoliubov averaging method known in the theory of ordinary differential equations 
(ODEs). As a matter of fact, the modulation equations (1.9) can also be derived via a 
rigorous multiple-scale perturbation procedure, an appropriate nonlinear analog of the 
classical WKB method involving an asymptotic expansion in a formal small parameter 
characterizing the scale separation. The modulation equations then naturally arise 
as the solvability conditions, in the space of periodic functions, for the first-order 
approximation. The leading order approximation yields the periodic solution itself. 
This approach, however, lacks directness and explicitness of the averaging method 
and can involve rather cumbersome calculations, especially in non-integrable cases. 

It should be stressed that the Whitham system (1.9) describes modulations of 
nonlinear periodic waves. In the hyperbolic case, the matrix A(u) generically has 
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N distinct real eigenvalues, which are the system’s characteristic speeds, sometimes 
referred to as nonlinear group speeds. In the limit of vanishing amplitude, two of 
these characteristic speeds merge together and degenerate into the linear group ve¬ 
locity uj' 0 (k), where lo = uJo(k) is the dispersion relation of the linearized original 
evolution equation. The splitting of the linear group velocity for modulated waves of 
finite amplitude is one of the most important results of the Whitham theory (see [98] 
Cli. 15.4). This effect leads to the expanding structure of DSWs. 

In the majority of dispersive hydrodynamic problems involving DSW formation, 
the modulation systems are hyperbolic, ensuring modulational stability of the DSW 
wavetrain. Hyperbolicity of the modulation system also enables use of the well de¬ 
veloped theory of quasilinear systems, such as the method of characteristics, which is 
not available for the original, dispersive evolution equation. 

By applying the Whitham method, the asymptotic description of a DSW ef¬ 
fectively reduces to the integration of the quasilinear modulation system (1.9) with 
certain initial/boundary conditions. These are formulated in terms of a continuous 
matching of the averaged, mean flow in the DSW region with the smooth external flow 
at boundaries that are determined along with the solution [40, 22], This nonlinear 
free boundary problem and its extensions are now commonly known as the Gurevich- 
Pitaevskii (GP) matching regularization. Importantly, the Whitham equations sub¬ 
ject to the GP matching conditions admit a global solution describing the modulations 
in an expanding DSW. The simplest, yet very important and representative, example 
of such a solution was obtained in the original paper [40] for a prototypical problem 
of the dispersive regularization of an initial step, the Riemann problem, for the KdV 
equation. The DSW modulation solution was found to be a rarefaction wave solution 
of the Whitham equations. For a more general case of monotone decreasing initial 
data, the existence of a global modulation solution describing the KdV DSW was 
proved in [93]. 

The GP asymptotic DSW construction procedure for the case of the Riemann 
problem has been extended to non-integrable systems by a priori assuming the exis¬ 
tence of a rarefaction wave solution to the Whitham modulation equations [22]. This 
assumption and the GP matching conditions yield the key, physical features of the 
DSW including the leading and trailing edge speeds, solitary wave edge amplitude, 
and small amplitude wave edge modulation wavenumber (see figure 1.1). 

Thus, within the context of the Whitham modulation approach, the dispersive 
regularization of wave breaking in a hydrodynamic system is implemented by the 
introduction of additional hyperbolic conservation laws (1.7), (1.8) augmented by 
certain free boundary matching conditions. This procedure is a dispersive counterpart 
to the classical regularization of viscous shocks, where the inviscid conservation laws 
are augmented by the Rankine-Hugoniot relations and entropy conditions. 

The modulation description of a DSW is formal in the sense that it is based on 
certain assumptions about the asymptotic structure of the DSW solution and thus, 
does not resolve the initial DSW development. Its rigorous validation requires the 
availability of the exact solution, and thus, is only possible for completely integrable 
systems. The modulation regularization procedure has been rigorously justified for 
fundamental integrable models such as the KdV and clefocusing Nonlinear Schrodinger 
equations in the framework of the zero dispersion limit of the Inverse Scattering Trans¬ 
form (1ST) developed by Lax, Levermore and Venakides [65, 94, 48]. In non-integrable 
cases, the validation of modulation theory is usually made via comparison with nu¬ 
merical solutions of the original dispersive initial or boundary value problem. We 
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note that, even for integrable equations, the direct application of modulation the¬ 
ory to the description of a DSW often has a number of advantages to the complex 
1ST based approaches. Being quasilinear hyperbolic equations, the modulation equa¬ 
tions do not have any global restrictions imposed on the solution, e.g., rapid decay 
as | a: | —> oo or spatial periodicity. The application of rigorous methods based on the 
1ST usually involve such limitations. Also, the applicability of the Whitham method 
to non-integrable equations makes the modulation regularization especially attractive 
for physical applications. 

Yet, the modulation approach has its own limitations. In particular, the results 
of modulation theory have an inherent long-time asymptotic character and do not 
yield important information about the initial stages of DSW development. A short 
time asymptotic theory has recently been developed by Dubrovin [19], which seeks to 
identify the universal oscillatory structure of dispersive hydrodynamic wave breaking. 
Another limitation to modulation theory is the fundamental assumption about the 
DSW structure as a locally periodic wave, which can be violated in some regions of 
the (x, t)-plane, in particular near the DSW edges [34, 35]. Other approaches such as 
rigorous Riemann-Hilbert steepest descent methods [16], [21] as well as the combina¬ 
tion of 1ST with matched asymptotic expansions [3] can yield additional information 
about the detailed structure of the oscillations. However, for the majority of applica¬ 
tions, the results of modulation theory usually provide the sought for qualitative and 
quantitative information of DSW behavior. The principal feature of the modulation 
approach is that it provides a unified framework for the analysis of dispersive regu¬ 
larization of nonlinear wave breaking singularities within the well established context 
of the theory of one-dimensional quasilinear systems. 

1.3. Overview of this work. In this paper, we compare solutions of the purely 
dispersive rnKdV equation to solutions of the diffusive-dispersive mKdVB equation 
with comparable initial data. We find quite remarkable similarities, but also striking 
differences in the structure of solutions. It is also of interest that the comparison 
depends on the sign of the dispersion coefficient /x; for both mKdV and mKdVB 
equations, the sign of p affects the structure of some solutions. In each case (p > 
0 , p < 0), we classify solutions of initial value problems with jump initial data: 


u(x, 0) 


U-, x < 0 
M+, X > 0 


( 1 . 12 ) 


according to the values of the constants u_, u+. We refer to these special initial value 
problems as Riemann problems. In our comparison of solutions of Riemann problems, 
we consider times long enough that sustained wave patterns have emerged. These are 
the long-time asymptotics available explicitly due to the simplicity of the mKdVB 
equation. 

In the purely dispersive case, solutions typically involve DSWs, expanding, slowly 
modulated wavetrains that provide for a dynamic transition between two different con¬ 
stant values of u. As described in the previous section, the asymptotic description of 
DSWs can be achieved by an averaging procedure that results in a system of quasi- 
linear first-order partial differential equations supplemented with initial/boundary 
conditions. DSWs resulting from Riemann problems are described by similarity so¬ 
lutions of the modulation equations. In contrast, solutions of the diffusive-dispersive 
equation involve traveling waves (TWs), solutions of an ordinary differential equation. 
The TWs approximate shock waves, which are discontinuous weak solutions of the 
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conservation law 


u t + ( u 3 ) x = 0. (1.13) 

This equation has characteristic speed 3it 2 , the speed of small disturbances, analogous 
to the sound speed in a compressible gas. Shock waves travel at the speed s = 
u 2 ,+Uj r U-+u 2 _ depending on the left and right limits, given by the Rankine-Hugoniot 
condition. Consequently, weak shocks (for which |w+ — u_| is small), travel at roughly 
characteristic speed. 

One of the most striking comparisons between the two equations is that for /i > 0, 
the mKdV and mKdVB equations exhibit special TWs termed undercompressive (in 
the language of [47]) because the corresponding shock wave 


u(x,t ) 


U-, X < st 

U+, X > st 


is subsonic both ahead of and behind the wave: 


(1.14) 


s<3u%. (1.15) 

The undercompressive TWs for mKdV have been known in the literature as kinks. 

Undercompressive shocks are non-classical in the sense that they do not satisfy 
the Lax entropy condition. They are called undercompressive by analogy with gas 
dynamics, in which classical or Lax shocks involve a jump in pressure associated with 
compression of the gas as a shock wave passes a fixed observer. When a shock fails 
to satisfy the Lax entropy condition, it is termed non-classical provided it has some 
relevance for the equation or application under consideration. For scalar equations, 
the only non-classical shocks are undercompressive in the sense of being either su¬ 
personic (see [91]) or subsonic on both sides of the wave. However, for systems of 
conservation laws, such as a system of magnetohydrodynamics [99], there can be a 
rich collection of non-classical shocks that have physical relevance, including both 
undercompressive and overcompressive shocks. See [88] for examples of 2 x 2 systems 
with both undercompressive and overcompressive shock waves. 

For /i < 0 there are no undercompressive shocks. In this case, the effect of 
the non-convex flux on the diffusive-dispersive dynamics (1.4) is the occurrence of a 
double-wave structure known as a shock-rarefaction. In the limit v —> 0, the shock 
wave in the shock-rarefaction represents a single-sided contact discontinuity propagat¬ 
ing with characteristic speed [15]. Contact discontinuities find their purely dispersive 
(v = 0) counterparts in the form of non-standard, contact DSWs whose modulations 
are described by special, double-characteristic expansion-fan solutions of the mKdV 
modulation equations. Such solutions become possible owing to nonstrict hyperbol- 
icity of the mKdV-Whitham modulation system. 

The shock wave theory for the mKdVB equation has been constructed in [47] by 
analyzing traveling wave solutions. Generalizations to different equations and to the 
Cauchy problem for equation (1.1) using wave front tracking are detailed in the book 
by LeFloch [68], 

There have been several recent developments in the modulation theory for mKdV 
DSWs. The two types of DSW solutions to the mKdV equation with p < 0 were 
constructed in [75] using the modulation equations obtained in [18]. These included 
the usual, KdV-type (convex flux) DSW solutions as well as non-standard sinusoidal 
undular bores, which we identify here as contact DSWs. The Riemann problem for the 
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closely related defocusing complex mKdV equation was considered in [61]. Elements 
of the DSW theory for the mKdV equation with n > 0 appear in [11] in the context 
of collisionless plasma physics. A complete classification of the Riemann problem 
solutions for the Gardner equation, a combined KdV-mKdV equation reducible to 
mKdV, was constructed in the recent paper [54] for both signs of ji. Application of 
the results for the Gardner equation to the resonant generation of internal undular 
bores in stratified fluid flows over topography was considered in [55]. 

Throughout this review, we incorporate direct numerical solutions in order to 
provide a visual rendering of the underlying mathematics. Here we fully resolve 
the small scale oscillatory behavior by utilizing a sufficiently small spatio-temporal 
discretization (see Appendix B). Note that numerical methods producing large scale, 
weak solutions of diffusive-dispersively regularized non-convex conservation laws that 
avoid the small scale discretization are an active area of research [70]. 

In §2, we execute a derivation of the KdVB and mKdVB equation from a general 
system of Eulerian equations for isentropic gas dynamics with non-convex pressure 
law regularized by a combination of viscous dissipation and a generic, third-order 
dispersion, using nonlinear multiple scale expansions. Similar derivations for partic¬ 
ular cases can be found in [50, 2, 96, 85]. The main purpose of our derivation is 
to identify the generic conditions under which the mKdVB model dynamics persist. 
The condition for the asymptotic mKdVB dynamics to dominate is shown to coin¬ 
cide with the condition for the loss of genuine nonlinearity of the Eulerian system 
in the non-diffusive, dispersionless limit. The important result of this section is the 
expression for the dispersion coefficient in the resulting mKdVB equation in terms 
of the pressure law and linear dispersion relation of the original non-convex Eulerian 
system. The sign of this coefficient is given by (2.13) and plays the definitive role in 
determining the qualitative structure of solutions of the Riemann problem for both 
mKdV and mKdVB equations. 

Section 3 is devoted to the comparison of classical Lax shocks with their purely 
dispersive counterparts by reviewing and comparing Riemann problem solutions of 
the KdVB and KdV equations. The main feature of a DSW, which is a key element 
of the KdV Riemann problem solution, is its expanding nature characterized by two 
distinct speeds, in sharp contrast with the steady structure of diffusive-dispersive 
shocks propagating with a single shock speed and described by TW solutions of the 
KdVB equation. The notion of a classical DSW is introduced via an analog of the Lax 
entropy condition formulated in terms of the characteristic speeds of the Whitham 
modulation system. It is shown that, for the KdV equation, the DSW admissibility 
condition coincides with the classical shock entropy condition, although the associated 
Rankine-Hugoniot condition is not applicable to the purely dispersive case due to the 
unsteady, expanding nature of the DSW. 

In §4, a unified description of the Riemann problem solutions for the mKdVB 
equation with positive and negative dispersion is presented following the results of 
[47]. The detailed structure of TW solutions approximating classical and non-classical 
(undercompressive) shocks is studied, distinguishing between oscillatory and non- 
oscillatory regimes. 

In §5, the full classification of Riemann problem solutions to the mKdV equa¬ 
tions with positive and negative dispersion terms is constructed in the framework 
of Whitham modulation theory, following the results of the previous study [54] per¬ 
formed for the Gardner equation. Families of classical and non-classical DSWs are 
then constructed and the admissibility conditions (5.20) for classical DSWs are in- 
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troduced. These conditions include the entropy-analog causality condition and the 
convexity condition implying a single-wave regularization. In contrast with the KdV 
case, the classical DSW admissibility conditions for the mKdV equation differ from 
the counterpart conditions obtained for rnKdVB classical shocks in the zero-diffusion 
limit. The distinct cases, where one or both of the classical DSW admissibility condi¬ 
tions fail, fully determine the boundaries in the U--U+ plane of initial data between 
regions of solutions with qualitatively different behaviors. Non-strict hyperbolicity of 
the modulation equations is shown to play the key role in the occurrence and structure 
of non-classical DSW solutions such as kinks and contact DSWs. 

Section 6 is central and provides a bridge between the two theories reviewed in §4 
and §5. It is devoted to a detailed comparison of the Riemann problem solutions to the 
rnKdVB equation with the solutions to the rnKdV equation with comparable initial 
data. This comparison reveals a number of non-trivial parallels and contrasts between 
the solutions of the two equations. We produce characteristic diagrams for diffusive- 
dispersive and dispersive regularizations. These diagrams, in particular, clarify the 
complex structure of non-classical rnKdV solutions and further elucidate the singular 
nature of the zero-diffusion limit. One of the manifestations of this singular nature is 
that the mapping between solutions of the rnKdVB and mKdV Riemann problems is 
not one-to-one. The ambiguities are resolved in the small-diffusion regime by identi¬ 
fying relevant time scales, including a transient window during which the dispersion- 
dominated dynamics evolve to accommodate the diffusive-dispersive balance. Finally, 
we show how the two contrasting analytical approaches to the description of diffusive- 
dispersive shocks and DSWs can be reconciled for the regime of small diffusion in the 
framework of a perturbed modulation theory. 

In §7, we draw conclusions from our study and outline future perspectives. The 
appendix contains a description of the numerical method used to compute solutions. 

2. KdVB and rnKdVB equations as weakly nonlinear approximations 
of diffusive-dispersive Eulerian hydrodynamics. The KdVB and rnKdVB equa¬ 
tions are prototypical equations, which often arise as unidirectional, weakly nonlinear 
long-wave approximations of more general diffusive-dispersive equations. For example, 
the Euler equations for isentropic gas dynamics in one space dimension, regularized 
by a combination of viscous dissipation and a generic third-order dispersion take the 
form, 

Pt + {pv) x = 0, (2.1) 

v t + vv x + p~ Y [P{p)]x = vv xx + [D(p, u)]a,. 

Here p > 0 is the density and v the velocity of the fluid, although these labels do not 
necessarily reflect the actual meaning of the corresponding physical entities (e.g. in 
nonlinear optics). In Eulerian dispersive fluids, the pressure law P(p) is an increasing 
function, P'(p) > 0. The parameter P > 0 is the viscosity and D(p,v) is a second 
order differential operator. We mention that the inviscid version of system (2.1) for 
certain classes of the operators [D(p,v)] x is sometimes called the Euler-Korteweg 
system [ 6 , 7]. 

The linear dispersion relation for system (2.1) has the form 

w(fc, po, v 0 ) =v 0 k-^k 2 ± u 0 (k, p 0 ), (2.2) 

obtained by linearizing ( 2 . 1 ) about the uniform background state p = po, v = vq 
and seeking the solution in the form of a traveling wave proportional to exp[i(kx — 
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cot)]. The ± signs in (2.2) correspond to the right- and left-propagating wave modes 
respectively. 

The key property of the dispersion relation we need is determined by the long¬ 
wave expansion of u>o(k, po), which we assume generically has the form 

cJo(k, p 0 ) = c 0 k + pk 3 + o(k 3 ), 0 < k <C 1, p ^ 0, (2.3) 

where Cq = c(po), c(p) = \JP'(p) being the long wave speed of sound. The sign 
sgn p = sgn{LL)g(/c,po)} is called the sign of dispersion and plays an important role in 
our consideration. The aim of this section is to outline the derivations of the KdVB 
and rnKdVB equations from the generic system (2.1) of isentropic Euler equations 
modified by a combination of viscous and dispersive corrections in order to identify 
the conditions under which the model dynamics, described in the subsequent sections, 
persist. 

We consider uni-directional, weakly nonlinear long-wave approximations of the 
diffusive-dispersive Euler equations (2.1) by introducing the following multiple scales 
expansions (see, e.g., [50, 2]). 


P — PO + £ Pl + P'2 + • • •, 

v = v 0 + £V\ + £ 2 v 2 + ■ • •, 


(2.4) 


where p t = Pi(x,r), v t = Vi(x,r),i > 1, and 

X = e p [x - Ut\, t = £ q t, 0 < £ < 1, p,q> 0. (2.5) 


To keep nonlinearity, dispersion and dissipation in balance we also need to require 
v = 0(e p ). To this end, along with (2.5) we introduce 


v = 


v 

2£P' 


( 2 . 6 ) 


To be definite we consider rightward-propagating waves (the plus sign in (2.2)) whose 
speed in the long-wave limit is U = vq + cq. We also assume, for i> 1, that pi , Vi —> 0 
as x ~:► +oo. 

1. KdVB approximation 

The choice p = 1/2, q = 3/2 yields, assuming the long-wave dispersion behavior 
(2.3), the KdVB equation for the first order corrections tq, p\, 

Mr + PM v i)x = "Mxx + P( v i)xxx . ( 2 - 7 ) 


where 


o = , Po p "(Po) 1 , Pqc'(po) 

1 2P'(po) + co 


( 2 . 8 ) 


Also, we get 


P o 

Pi = —vi . 
co 


2. mKdVB approximation 
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The KdVB approximation (2.7) is valid as long as |/3| = 0(1), e -A 0 in (2.8), 
which is always the case for systems with convex pressure law, P"(p) > 0. However, 
if P"{p) < 0 for some interval of p, then there is the possibility that the nonlinear 
coefficient f3 could be zero. This occurs if the background density po satisfies the 
criticality condition 

(p 2p '(p)Y\ p = P0 = 0, or, equivalently (pc(p))'\ p=po = 0. (2.9) 

On the other hand, if /3 remains non-zero, then the hyperbolic limit of equations (2.1) 
are genuinely nonlinear in the sense of Lax [64] (see, e.g., [15]). Genuine nonlinearity 
is most succinctly formulated as convexity of the pressure expressed as a function 
of the specific volume w = p _1 : V"{w) ^ 0 where V{w) = P(l/w;). Our main 
interest will be in non-convex systems , for which the genuine nonlinearity condition 
is violated at exactly one point. In that case, third-order terms become necessary in 
the expansions (2.4) to account for the higher-order nonlinear term in the resulting 
equation for v\. Then balance is achieved with new scaling parameters p = 1 and 
q = 3 in the transformation (2.5) (see, e.g., [96], [85]). Substituting the new multiple 
scales expansions (2.4), (2.5) into (2.1) and equating like powers of e up to e 4 , we 
obtain the mKdVB equation for the first order corrections, 

(wi)r + 7 (Ti) 2 (Ti)x = v(vi)xx + K v i)xxx - Pi = ~ v 1- 

Co 


Here, 


and 


Introducing 


X = z[x — (vq + c 0 )t], r = £ 3 t, 


7 = ^3 k '"(Po) . wh ere «(p) = pP(p). 


u — 



11 = /isgn 7 , 


( 2 . 10 ) 


( 2 . 11 ) 


( 2 . 12 ) 


and replacing r by t, \ by sgn(q)a:, we arrive at the standard form (1.4) of the mKdVB 
equation. 

If /3 in (2.8) is not strictly zero but O(e), the slightly modified long-wave multiple 
scales expansions (2.5) with p = 1, q = 3 lead to the mixed, KdV-mKdV-Burgers 
equation, with both quadratic and cubic nonlinearities present. The zero-diffusion 
version of this equation is often called the Gardner equation. See [85] and [53] for the 
multiple scales derivations of the Gardner equation for ion-acoustic waves in collision- 
less plasmas with negative ions and nonlinear polarization waves in two-component 
Bose-Einstein condensates respectively. 

As mentioned above, the structure of some solutions to the Riemann problem for 
the mKdV and mKdVB equations depends crucially on the sign of the dispersion co¬ 
efficient p, which is defined in terms of the key nonlinear and dispersive characteristics 
of the original Eulerian system (2.1) as 


sgnp = sgnln'”(po)u k k(k, po)} ■ 


(2.13) 
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In systems with non-convex dispersion law u)o(k,po), in addition to non-convexity 
of the pressure law, there is the possibility of a sign change in p due to a range 
of values of k involved in the Riemann problem solutions, which can lead to further 
complexity in their structure. See [72] for the Riemann problem analysis of a dispersive 
hydrodynamic equation with convex hydrodynamic flux but non-convex dispersion 
law. However, the possible change of dispersion sign is likely to occur for sufficiently 
short waves and thus, may not affect the long-wave approximation. 

3. Lax shocks and their dispersive counterparts. The shock wave (1.14) 
is a weak solution of the conservation law (1.1) if it satisfies (1.5), from which the 
Rankine-Hugoniot condition 

-s(u+-u-) + f(u+)-f(u-)= 0 (3.1) 

can be deduced. Such a solution is a Lax shock if it satisfies the Lax entropy condition 
relating the characteristic speeds f'[u±) to the shock speed s = (f(u+) — f(u-))/(u+ — 
u -) : 


/'(«+) < s < (3.2) 

For the inviscid Burgers equation 

ut + \{u 2 ) x = 0, (3.3) 

Lax shocks are characterized by s = ^(u++U-), with u+ < it_. In the next subsection, 
we show that shock wave solutions of (3.3) are Lax shocks if and only if there are 
traveling wave solutions of the KdVB equation 

Ot + 2 ( u )x = VU X x "f P^xxx (3.4) 


with v > 0. 

In §3.2, we describe the purely dispersive counterpart of the Lax shock for v = 0 in 

(3.4) , the KdV equation. The main point we wish to emphasize here is that solutions 
of the KdV equation with initial data a Lax shock are quite different in character 
from those of the KdVB equation, and possess a structure quite distinct from the 
Rankine-Hugoniot condition or the traveling wave solutions of the KdVB equation 

(3.4) with v > 0. 

3.1. Traveling wave solutions approximating shocks. We wish to consider 

(3.4) with p ^ 0. Since the change of variables x — > —x, u —> — u leaves the equation 
unchanged if p is replaced by —p, we can take p < 0 without loss of generality. Then 
if we seek TW solutions in the form u(x, t) = u((x — st )/ vTa'T) with u(±oo) = u± , we 
see that such solutions approach the shock wave (1.14) as \p\, v —► 0 + and it remains 
only to demonstrate the existence of the TW solutions. 

We find after integrating once that TW solutions satisfy a second order equation 


—s(u — U-)+^(u 2 — u 2 _) = Su l ~u", d = v/\f\p\, (3-5) 

where £ = (x — st)/^/\p\, ' = d/d^, and we have dropped the tilde from u. If u(£) 
is a solution for a fixed S, then as p —> 0, we have v —> 0 also. Hence, the traveling 
wave approaches a shock wave, which is in fact a weak solution of the conservation 
law (3.3), i.e., (3.4) with p = v = 0. Specifically, assuming that v! and u" approach 
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u 

A 


<5 = " 
0 < 
vf 



u 


Fig. 3.1. Phase plane depicting KdVB TW solutions connecting u_|_ < U— when // < 0. The 
dashed curve is a monotone TW with strong diffusion. The solid, spiraling curve is an oscillatory 
TW with weak diffusion. Both are contained inside the homoclinic orbit of the zero diffusion KdV 
solitary wave solution. The filled region is the union of all KdV periodic orbits spanned across a 


DSW. 


zero at ±oo, the Rankine Hugoniot condition is satisfied: s = \{u + +it_). The second 
order equation (3.5) can be written as a first order system 


u' = v, v' = Sv — u 2 — u 2 _) + s(u — U-). 


For S > 0, the phase portrait of this system has two equilibria (u, v) = (u±, 0), one of 
which is a stable saddle point and the other is an unstable node or spiral, depending 
on whether the eigenvalues of the linearized system 



are real (<5 > V2A) or complex (<5 < \/2A), respectively, where A = |u_ — u+\ is the 
jump across the TW (see figure 3.1). An analysis of the phase portrait reveals that 
the saddle point is connected to the second unstable equilibrium by a heteroclinic 
orbit, and this corresponds to the TW we seek if and only if u + < u_, namely the 
Lax entropy condition. (Note that for /i > 0, the previously unstable node or spiral is 
stable). The structure of the TW changes from a monotonic to oscillatory profile as S 
is decreased across the critical value v / 2A- Further analysis shows that for 0 < <5 <C 1, 
the leading amplitude of the oscillatory profile is approximately a = |A [49]. For 
<5 = 0, however, the unstable equilibrium becomes a center, i.e., with imaginary 
eigenvalues, and a TW connecting u + and u_ no longer exists. Instead, there is a 
homoclinic orbit connecting the saddle point to itself, and the corresponding TW is 
a KdV solitary wave. Traveling wave solutions to the KdVB equation were analyzed 
in [33, 49, 9], 

When u + > w_, the inviscid Burgers equation (3.3) has a rarefaction wave solution 



—, U-t < x < u+t 

U- x > u+t. 


(3.6) 
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U +A 


U + = U- 


u + A 


u + = u 


-► 

u _ 


a) p < 0, 0 < v < y/2\fi\ 


b) p > 0, 0 < v < y/2/j, 


Fig. 3.2. Classification of KdVB Riemann problem solutions, (a) Negative dispersion; (b) 
Positive dispersion. 

This solution requires that u + > . Although rarefaction waves have corners at the 

edges x = u±t, these irregularities are smoothed by the dissipative and dispersive 
terms, so that apart from this small correction, the rarefaction waves appear for the 
full KdVB equation. 

We can now classify the long-time behavior of solutions of the Riemann problem 
for (3.4). For u- > u + , the solution approaches (as t -A oo) the traveling wave 
connecting u- to u+, whereas for < u+, the solution is approximately a rarefaction 
wave. This classification is depicted in figure 3.2. 

3.2. KdV equation and DSWs. As remarked at the end of §3.1, when v = 0 
in (3.4), there is no traveling wave corresponding to a heteroclinic orbit in the phase 
portrait, hence the Rankine-Hugoniot condition for a shock wave is not satisfied. In¬ 
stead, there is a homoclinic orbit connecting the saddle point to itself (a soliton), and 
there are periodic orbits representing periodic traveling wave solutions of (1.3). In 
particular, the approximation of a shock by traveling waves is invalid when v = 0. 
This could lead one to the conclusion that a certain amount of dissipation is nec¬ 
essary to sustain a shock wave transition. Indeed, the classical theory of undular 
bores by Benjamin and Lighthill [5] and the theory of collisionless shocks in plasma 
due to Sagdeev [86, 87] and others (see, e.g., [33]) are based on dissipative-dispersive 
models, essentially reducible to the KdVB equation (3.4). It was not until Gurevich 
and Pitaevskii’s pioneering work [40] that it was realized that there is a qualitatively 
different type of shock wave that develops in the total absence of dissipation. These 
conservative shocks are called dispersive shock waves (DSWs) and represent expand¬ 
ing, modulated nonlinear wavetrains connecting the constant states u_ and u.y. 

The dispersion relation for the KdV equation (1.3) on the background uq is 
ui(k,uo) = uok + /ifc 3 . It will be convenient to first consider the case of negative 
dispersion, (i < 0. DSW solutions of the KdV equation are asymptotically described 
by periodic solutions (cnoidal waves) 



(3.7) 
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where rq < ri < r 3 are parameters, which remain approximately constant on the 
wavelength/period scale Ax, At ~ a/H; dn(£,m) is a Jacobi elliptic function, 


x — Ut + Xq 

V\iA 


u = 


1 

3 


(f 1 + r 2 + r 3 ), 


i ’2 - n 

m =-. 

r3 r 1 


(3.8) 


Here Xo/^/\ji\ is the initial phase and the modulus m satisfies 0 < to < 1. As to is 
varied from 0 to 1, the corresponding periodic orbits sweep out a phase plane region 
enclosed by the solitary wave homoclinic orbit when m —> 1~ (see figure 3.1). The 
limit to -A 0+ corresponds to the linear, harmonic wave limit. Slow, spatiotemporal 
modulations r,(a;,t) of the function (3.7) occur on a scale substantially larger than 
\J |p | and provide the required transition between the two distinct constants and 
u + so that to = 0 at one edge of the DSW and to = 1 at the other edge. Thus the 
DSW completely fills the phase plane region enclosed by the solitary wave orbit. In 
figure 3.1, we see that the behavior of the KdVB TW for weak diffusion 0 < S < \/2 
spirals through, but does not fill, this same region. 

The modulations are described by a system of three quasilinear, first order (hy¬ 
drodynamic type) equations derived originally by Whitlram [97] by averaging three of 
the conservation laws of (1.3) over the periodic family (3.7) resulting in 

dv ■ dv' 

-gij- + V i (r 1 ,r 2 ,r 3 )-^ = 0, * = 1,2,3, (3.9) 

i.e., r 1 , v 2 , r 3 are Riemann invariants for the modulation equations. The characteristic 
speeds V\ < V 2 < V 3 are found as (see (1.11), (3.8)) 

Vi = u ~vk' ( 3 - 10 ) 

6 dri 

where 

L = 2K(m)J —-—, (3.11) 

V r 3 - r 1 


K(to) being the complete elliptic integral of the first kind. The quantity LsJ |7*[ is the 
wavelength of the periodic wave (3.7). The solution r (x,t) of the Whitlram system 
(3.9) with appropriate initial or boundary conditions is inserted into the cnoidal wave 
(3.7) to obtain the asymptotic solution for a DSW. We note that the initial phase 
Xq in (3.8) remains undetermined within Whitlram theory so the outlined solution 
describes the local structure of a modulated wavetrain up to a spatial shift ~ \/|/i 
within the traveling wave period. 

We now present the explicit expression for , r' 2 , 7 * 3 ), which plays a particularly 
important role in the DSW theory. 


V 2 (r 1 ,r 2 ,r 3 ) 


g(»T + r 2 + r 3 ) 


2 (1 — to)K(to) 

3 (r2 "" Tl ’ E(to) - (1 - to)K(to) ’ 


where E(to) is the complete elliptic integral of the second kind. Of particular relevance 
are two limits of V 2 corresponding to the harmonic and soliton limits of the traveling 
wave (3.7). In the harmonic limit, the modulus to = 0 (i.e. r 2 = rq, see (3.8)) and 
the characteristic speeds V 2 and V\ merge together 


to = 0 : V 2 (r 3 , ri,r 3 ) = Hi(ri, rq, r 3 ) = -r 3 + 2r 3 = VL . 


(3.12) 
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One can also show that in this limit the velocity V 3 (?q, rq, rq) = r 3 . Thus, in the 
harmonic limit the Whitham system (3.9) reduces to a system of two equations, one 
of which is the inviscid Burgers (dispersionless KdV) equation (3.3) for r 3 . 

In the opposite, soliton limit m = 1 (i.e., r 2 = r 3 ), we have a similar degeneracy 
but now the merged characteristic speeds are V 2 and V 3 

m = 1 : V 2 (n,r 3 ,r 3 ) = V 3 (n,r 3 ,r 3 ) = ^(2r 3 + rq) = V+ , (3.13) 

while the remaining velocity Vi(rq, rq, rq) = rq yields the inviscid Burgers equation 
for rq. 

The described reductions of the KdV-Whitham system (3.9) in the harmonic and 
soliton limits enable one to construct the solution r (x,t) = (ri(x,t),r 2 (x,t),r 3 (x,t)) 
that provides a continuous matching of the period-mean field, 

u(ri,r 2 ,r 3 ) = - / u(£)d£ = n + r 2 - r 3 + 2(r 3 - rq) —, (3.14) 

L J 0 K (to) 

with the distinct constant states u_,u+. This results in u(rq, 7 q,rq) = r 3 = U- at 
some x = X-(t) (the DSW’s trailing, leftmost edge) and «(rq,rq,rq) = rq = u + at 
some x = x+(t) (the leading, rightmost edge). The positions of the DSW edges are 
free boundaries defined by the merged characteristics dx±/dt = V± evaluated for the 
relevant modulation solution. The outlined matching regularization problem is due 
to Gurevich and Pitaevskii [40], and its solution r (x,t) describes slow modulations 
of the periodic solution (3.7) in a DSW. We note that, while the modulation theory 
defines the “carrier” cnoidal wave solution (3.7) only up to an arbitrary phase shift 
Ax ~ macroscopic DSW parameters (the amplitude, the wavenumber, the 

edge speeds, etc.) are fully determined by the solution of the Gurevich-Pitaevskii 
problem. The computation of the exact phase Xo/\f\jj\ in the DSW solution would 
generally require invoking a higher-order approximation, beyond the original leading- 
order Whitham theory. Despite this, for some classes of initial data the phase can 
be explicitly expressed in terms of the Whitham modulation parameters rq [16], [34], 
[24]. For the case of the Riemann problem considered in this paper, the DSW phase 
was derived in [3]. 

The KdV-Whitham system (3.9) is strictly hyperbolic and genuinely nonlinear 
[71], i.e. for all r£l 3 one has: (i) 7 q 7 ^ rq <=$■ V) 7 ^ 14 (strict liyperbolicity), and 
(ii) dVj/drj ^ 0, j = 1,2,3 (genuine nonlinearity). Indeed, the merged characteristic 
speeds (3.12) and (3.13) in the harmonic and soliton limits respectively define a regular 
characteristic since the order of the modulation system reduces to two in both limits, 
and thus strict hyperbolicity and genuine nonlinearity of the modulation system are 
preserved in these limits. As a result, the dispersive regularization of the initial step 
(1.12) is described by a self-similar simple-wave modulation of (3.7) which is the 
rarefaction wave solution of the Whitham equations , 

X 

P < 0 : 7*1 =u + , r 3 = U-, V 2 {u + ,r 2 ,u-) = (3.15) 

Being a rarefaction fan (see figure 3.3c) solution (3.15) exists for all t > 0. It is 
defined in s-t < x < s+t, where s_ and s+ are the speeds of the trailing and the 
leading DSW edges respectively (see figure 3.3). They are found from (3.15) by 
setting r 2 = rq (m = 0 trailing, harmonic wave edge) and r 2 = r 3 (m = 1 
leading, soliton edge) to obtain (see (3.12), (3.13)) s_ = V 2 (u + ,u + ,u-) = u + — A 
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and s+ = V 2 (u+,u-,u-) = u+ + 1 A, where A = u- — u + is the jump across the DSW 
[40]. It is instructive to note that the speed of the harmonic wave edge coincides with 
the linear group velocity cok evaluated at some k > 0 while the soliton edge speed is 
determined by the phase velocity in the limit A; -A 0 [22]. The wave number k = 2ix/L 
(see (3.11)) in the DSW is a monotone function of x, which decreases as the DSW is 
traversed from the harmonic to solitary wave edge. For fj, < 0 one has dk/dx < 0. 

Lax and Levermore [65] showed that this solution of the Whitham equations 
describes the weak limit solution u(x,t,n) of the KdV equation (1.3) when /x —>■ 0, 
the weak limit being the period-mean (3.14) with n, r 2 , r 3 defined by (3.15). We stress 
that this weak limit does not coincide with the weak solution of the Riemann problem 
for the conservation law (3.3), which is the Lax shock (1.14) with s = («_ + u+)/2. 
This is due to the oscillations in the DSW having wavelength with scale ~ \/\lA so the 
limit as /x — > 0 is a singular one. Representations of the Riemann problem solutions 
in the x-t plane for the KdV equation with /i = 0 (the inviscid Burgers equation) and 
for fx —> 0 (the Whitham equations) are shown in figure 3.3. 



Fig. 3.3. Characteristics for the KdV Riemann problem solutions. (a,b) The inviscid Burgers 
equation p = 0. (a) U— > u+, classical shock; (b) U— < u+, rarefaction wave; (c-f) U— > u+, 
— >• 0, rarefaction wave of the Whitham equations describing the DSW modulation, (c) The 2- 
characteristic family T 2 of the Whitham equations exhibiting the rarefaction fan. (The characteristic 
corresponding to the classical shock is shown in red.) (d) The 1-characteristics, labeled Ti, approach 
x = s_t as t —»■ 00 (2-characteristics are shown as dashed lines for comparison), (e) The 3- 
characteristics, labeled T 3 , approach x = s+t as t —»■ 00 . (}) The Whitham Riemann invariants 
rj,j = 1,2,3 at a specific time. 


The general DSW admissibility conditions for a scalar dispersive equation with 
convex hyperbolic flux are [22], [24]: 

S- < s + > S- < s+ . (3.16) 

These conditions are also known as the DSW causality conditions. Similar to the clas¬ 
sical Lax entropy conditions (3.2), inequalities (3.16) ensure that the characteristics 




















20 


G. A. EL, M. A. HOEFER, AND M. SHEARER 


transfer initial data into the DSW region. The generalization of the DSW admis¬ 
sibility conditions (3.16) to bidirectional Eulerian dispersive hydrodynamics can be 
found in [22], [44]. We shall refer to DSWs satisfying the Lax-type conditions (3.16) 
as classical DSWs. 

For the KdV equation (1.3), the dispersionless (/i = 0) characteristic speeds are 
/'(«_) = U-, /'(«+) = u + , while speeds of the leading and trailing DSW edges 
(for ji < 0) are s + = u + + |A and s_ = u + — A. Thus, in the KdV case, the 
DSW admissibility conditions (3.16) are equivalent to the entropy condition u + < u_ 
for Lax shocks. However, we stress again that the Rankine-Hugoniot condition s = 
( u - + u +)/2 — u + + \ A is not applicable to a DSW. Moreover, the Lax shock speed 
satisfies s_ < s < s + (see the red line in figure 3.3(c)). 

As noted for the KdVB equation earlier, the transformation x —> —x, u —► —u, 
[i -A —fi leaves the KdV equation unchanged. Therefore, a reflection of the spatial 
and amplitude axes of the KdV DSW with negative dispersion, p < 0, is sufficient 
to describe the DSW with positive dispersion. Explicitly, the traveling wave solution 
and the modulation equations for the case of positive dispersion are obtained from 
(3.7) - (3.11) by applying the transformations 

[X 0 ^ fX S’ 0 U —y U, Vi ? —7*4—,;, V —^ V 4 _ 2 , i = 1,2,3. (3.17) 

We note that under the transformations (3.17), the modulus of the elliptic function 
transforms as m —> 1 — m. 

The change p < 0 —> fx > 0 alters the leading and trailing edge DSW speeds 
(s+ = V 2 = U+ + 2A, S- = V 2 (u+,u+,u-) = u+ + ^A when p > 0) but, 
perhaps more strikingly, the DSW structure is changed. The local description of the 
DSW according to (3.7) for /x < 0 takes the form of the KdV solitary wave 

u = u + + 2Asech 2 , 

as the parameter r 2 —> r 3 (m -A 1). Therefore, the leading edge of the DSW exhibits a 
positive, elevation wave with amplitude a = 2A. Whereas for positive dispersion /i > 
0, the trailing edge of the DSW exhibits a negative, depression solitary wave. Figure 
3.4 displays the two cases. The sign of the solitary wave, elevation or depression, 
defines the polarity p of the DSW. We denote a DSW that contains an elevation 
(p = 1) or depression (p = —1) solitary wave as DSW + or DSW - , respectively. 

The position of the solitary wave, at either the leading edge or the trailing edge, 
defines the orientation d of the DSW. We say that d = 1 if the solitary wave is at 
the leading edge, otherwise d = — 1. The orientation can be found by evaluating 
d = — sga[dk/dx] at the harmonic edge. This can be understood in very general 
terms as follows [38]. As already mentioned, the modulation wavenumber k decreases 
as the DSW is traversed from the harmonic to solitary wave edge. The orientation of 
the DSW is determined by a well-ordering of the DSW in the vicinity of the harmonic 
edge. The harmonic edge is at the left (trailing edge) if the group velocity increases 
with decreasing wavenumber, i.e., LOkk < 0 (p < 0). Conversely, the harmonic edge is 
at the right (leading edge) if ujkk > 0 (p, > 0). In both cases d = —sgnp i.e. the DSW 
orientation is uniquely defined by the sign of the dispersion. 

Thus for the KdV equation, both the DSW polarity and its orientation are de¬ 
termined by the sign of dispersion alone. As we shall see, for the mKdV equation, a 
change in the DSW polarity (but not the orientation) becomes possible without the 
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change of sign of the dispersion. These properties are inherited from the correspond¬ 
ing KdVB traveling wave solutions as indicated by a comparison between figures 3.2 
and 3.4. 

If U- < u+, the leading order solution of KdV with |/t] -C 1 represents a classical 
rarefaction wave (3.6). In contrast to the Riemann problem solution of the KdVB 
equation (3.4) with u_ < u+, one corner (left if < 0, right if n > 0) of the rarefaction 
wave is smoothed by linear dispersive oscillations whose wavelength ~ y/HA an d 
amplitude decays like f -1 / 2 , while the other corner is exponentially resolved to a 
constant [67]. The Riemann problem classification for the KdV equation is shown 
schematically in figure 3.4. 




a) n < 0, v = 0 


b) fi > 0, v = 0 


Fig. 3.4. Classification of the KdV Riemann problem solutions, (a) Negative dispersion; (b) 
Positive dispersion. 


Concluding this section, we summarize the contrasts between KdV DSWs and 
KdVB traveling shock wave solutions (see figure 3.5) with 0 < v <C y/2A/j,. 

1. The most striking contrast between the two structures is that the KdV DSW 
is an unsteady, expanding structure, and is characterized by two speeds which 
are defined as the speeds of the leading and trailing edge characteristics of the 
Whitham modulation equations. The KdVB TW describing the diffusively 
modified shock is characterized by a single classical shock speed determined 
by the Rankine-Hugoniot condition. 

2. The amplitude of the leading soliton in the KdV DSW is 2A while its coun¬ 
terpart in the shock solution of the KdVB equation with small v > 0 has 
amplitude 3A/2. 

4. mKdVB and non-classical shocks. In this section, we restrict attention 
to the diffusive-dispersive mKdVB equation (1.4) with v > 0. Our primary interest is 
to relate traveling waves to shock waves (1.14) satisfying the conservation law 

ut + ( u 3 ) x = 0. (4.1) 

From (3.2) we have that a shock wave solution (1.14) of (4.1) satisfies the Lax entropy 
condition if 


3u+ < s < 3u 2 . 


(4.2) 
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Fig. 3.5. Numerical solutions of the Riemann problem for KdV (dashed, green) and KdVB in 
the oscillatory (dash-dotted , blue) and monotonic (solid, red) TW regimes with the same step-like 
initial conditions and p = — 1. 


Here, 


s = u 2 + + u+u- + u 2 _ (4-3) 

is the wave speed. Thus, (1.14) is a Lax shock (sometimes termed a classical or 
compressive shock ) if and only if 

—sgn(u_)-u_ < sgn(u_)it+ < sgn(u_)u_ . (4.4) 

Note the asymmetry in the inequalities resulting from the combination of the Lax 
entropy condition and the Rankine-Hugoniot relation. 

When it + = —u_/2, the shock wave (1.14) satisfies the Lax entropy condition in 
a limiting sense, since s = 3u^_. That is, the shock speed is characteristic on the right. 
Such a shock is sometimes called a one-sided contact discontinuity. 

4.1. Traveling wave solutions. To characterize traveling wave solutions u = 
u(x — st) of (1.4) with speed s, we impose the far field boundary conditions 

u(±oo) = u±, u(± oo) = 0, ^(± 00 ) = 0. (4.5) 

Then it satisfies the ODE (dropping the tildes): 

— s(u — U- ) + u 3 — u 3 _ = vu' + \xv!'. (4.6) 

In particular, the boundary condition (4.5) at +00 yields the Rankine-Hugoniot con¬ 
dition 


— s(u + — U -) + — u 3 _ = 0, (4.7) 

giving the wave speed s in (4.3). 

The following definition depends on fixing the parameter S = zz/VT/rj anc ^ the s ig n 
of pt. A shock wave (1-14) is termed admissible if there is a traveling wave solution 
of (1.4) satisfying the far-held conditions (4.5). Consequently, the shock wave is 





































DISPERSIVE AND DIFFUSIVE-DISPERSIVE SHOCK WAVES 


23 


admissible if and only if the parameters s, ti+,ti_ are related by (4.3) and there is a 
solution of the ODE (4.6) satisfying ti(±oo) = u±. We note that the ODE is unchanged 
under the transformation u —» —u, U- —U-. Consequently, if (1.14) is an admissible 
shock with speed s, then so is the shock with u± replaced by —ti±, with the same 
speed s. 

Just as for the KdVB equation, we have a second order ODE that describes 
traveling waves ti(£),£ = (x — st)/\f\fi\, except that for the MKdVB equation, the 
sign of n is significant: 

sgn (n)u" + 5u' = u 3 — u ^ — s(u — ti_), 8 = v/ a/M- 

The corresponding first order system is 

v! = v, sgn(/i)t/ = —Sv + ( u 3 — it_) — s(u — w_). 

For S > 0, the phase portrait of this system has between one and three equilibria 
including (u,v) = (tt±, 0), where u+ satisfies (4.7). We are interested only in cases 
for which there are three equilibria, two of which may be a double root of (4.7). The 
nature of the equilibria depends on the eigenvalues 

A±(ti) = \ (—(sgn n )8 ± \/S 2 + (sgn/i)4(3ti 2 — s)^ , s = u 2 + uu- + v?_. (4.8) 

Consider the case of three equilibria with u = u± , uq satisfying u+ < u o < u _. Then 
shocks from u_ or u + (with x < st) to Uq (with x > st) satisfy the Lax entropy 
condition, since 3«g < s < 3 u±. 

In terms of the equilibria, the eigenvalues (4.8) imply that if p < 0, then (u, v) = 
(u±, 0) are unstable nodes or spiral nodes (if 8 2 +4(3itg — s) > 0 or S 2 + 4(3ug — s) < 0, 
respectively), and (ug,0) is a saddle point. By tracing the stable manifold of (u 0 , 0) 
backwards in time, we can establish that there is a trajectory from each of (u±. 0 ) to 
(uq. 0). Consequently, for /u < 0, all Lax shocks are admissible and there are no other 
admissible shocks. An example phase portrait is shown in figure 4.1(a). 

On the other hand, if p. > 0, then (u,v) = (u±, 0) are saddle points, and (uq, 0) 
is a stable equilibrium, either a node or a spiral. The saddle points at (u±, 0) may or 
may not be connected to (tig, 0), and can be connected to each other. See figure 4.1(b) 
for a phase portrait when both u± are connected to tio, corresponding to admissible 
Lax shocks. Figure 4.1(c) is a phase portrait showing a connection for a Lax shock 
but also an orbit connecting w_ to u + . In this case, there is no connection from ti_ 
to tig, even though these constants define a Lax shock. It is proved in [47] that an 
orbit joining saddle points u± exists if and only if it lies on an invariant parabola 
v = ±^(ti — ti_) (ti — ti_|_); an explicit calculation shows that this occurs precisely for 


|ti_| > 



ti+ = — ti_ + sgn(ti_) 



(4.9) 


for which s = u 2 _ — ^<5|ti_| + 1 8 2 , and tig = — sgn(ti_)^<5. These values help deter¬ 
mine which Lax shocks are admissible for /t, > 0, and which are not. The dynamics 
on the invariant parabola correspond to the traveling wave 

= ^ (u-~u+), £ = (x-st)/Vfi. ( 4 - 10 ) 


tt(^) = ^{u_+ti+ — |ti_—ti_|_ | tanh(A^)}, A 
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v v v 



Fig. 4.1. Representative phase portraits for traveling waves of the mKdVB equation. The 
solid curves correspond to Lax shocks connecting u± to uq. The dashed curve in (c) represents an 
undercompressive shock connecting u— to u+. 


To summarize, for \u-\ > a shock wave (1.14) from U- to u+ is admissible 

if either (i) u+ is between —sgn (u-) y ^-5 and U- or (ii) u+ = — it_ + sgn(M_)^d. 
In case (i), the shock satisfies the Lax entropy condition 3 u 2 + < s < 3 uf, and in 
case (ii) the shock is undercompressive, in the sense that it has the property (1-15) 
discussed in the introduction. In case (ii), other Lax shocks (for which rt + lies between 
— \u- and —sgn (u-)^-S) are not admissible. In the case |rt_| < ^r^<5, there is no 
undercompressive shock (for any u + ) and all Lax shocks (for which u+ lies between 
U- and — \u-) are admissible. 

When u + = — U-/2 , the shock wave (1.14) is a one-sided contact discontinuity. 
It is admissible if either fj, < 0, or n > 0 and |n_| < ^^6, the latter inequality arising 
from the existence of the saddle-saddle connection of Figure 4.1(c). 

4.2. The Riemann problem. To complete the description of solutions of the 
Riemann problem, we need to utilize not only constant solutions and admissible 
shocks, but rarefaction waves. These are continuous self-similar solutions of (4.1) 
that join constants u±, and are given by the explicit formula 

U-, X < ?>u 2 _t 

Suft < X < 3u+t (4-11) 

u+, x > 3 u\t. 

Rarefaction waves therefore propagate to the right, and u(x,t) has a single sign 
throughout the wave. Consequently, |u+| > |u_| > 0. Note that the graph of 
u = u(x,t ) as a function of x for fixed t > 0 is continuous, and smooth apart from 
corners where the parabola meets the constants u = u±. Using this information about 
shock waves and rarefaction waves, we can construct solutions of the Riemann problem 
for any u + ,u_ using only constant solutions, admissible shock waves (both classical 
and undercompressive), and rarefaction waves. The classification diagram for these 
solutions is presented in figure 4.2. The construction labeled R, with j = ±,0, 
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a) n < 0, S = vj vTH > 0 



Fig. 4.2. Solutions of the Riemann problem for the mKdVB equation with 5 = u/W\jj\ > 0. 
Line (T): = —u~/ 2. Line ©: u+ = —u~ + y/26/3. 


indicates a shock-rarefaction. This is a composite wave in which a shock from U- 
to —U-/2 has speed s = 2>u 2 _/2 that is characteristic on the right, at it = —U-/2. 
Consequently, to reach u+, beyond —U-/2 (meaning u+ < —«_/ 2 when u- > 0 for 
example), a faster rarefaction wave can be attached to the shock wave, yielding the 
single wave structure in which the rarefaction wave has a discontinuity at it’s trailing 
edge. The discontinuity is sometimes called a one-sided contact. 

In figure 4.2(a), in which /i < 0, we use the fact that there are no undercompressive 
shocks and all Lax shocks are admissible. In figure 4.2(b), we show the structure of 
solutions when /i > 0. The classification is more complicated, in that there are new 
regions in the diagram, in which the Riemann problem solution has two waves traveling 
with different speeds, separated by a widening interval in which u is constant. This 
structure is indicated in the figure with a vertical line between the wave letters. One 
wave is an undercompressive wave, indicated with the letter E, the other being a 
shock S or rarefaction R. Note that figures 4.2(a,b) inherit the symmetry u+ —> — u+, 
U- —> — U- from the governing mKdVB equation. 

Strictly speaking, figure 4.2 shows solutions of the Riemann problem for the con¬ 
servation law (1.13), using information about traveling waves with 0 < d < oo to 
distinguish those shocks that are admissible. But the conservation law can be re¬ 
garded as the expressing the leading order long-time behavior of solutions of the 
Riemann problem for the MKdVB equation as /i and v approach zero with <5 > 0 
held constant. In figure 6.6, we show numerical solutions of the Riemann problem in 
representative cases of each of the regions shown in figure 4.2. In that figure, shocks 
are replaced by smooth traveling waves, and rarefaction waves, though not precisely 
self-similar, are clearly visible. 

It remains to explain the superscripts ±,0 associated with shock waves. These 
indicate whether there may be oscillations (±) or not (0). As explained below in the 
context of purely dispersive shocks, S + indicates that the traveling wave terminates 
exponentially at an equilibrium, with an adjacent maximum. Similarly, S~ indicates 
that the largest oscillation involves a minimum before the traveling wave terminates 
exponentially at the equilibrium. 
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To explain these structures, we refer to the ODE (4.6). We can write the equation 
as a first order autonomous system 

v! = v, fiv' = —vv + f(u;u^.,s), f(u;u-,s)=u 3 —u 3 _—s(u — u-). (4.12) 

See figure 4.1 for representative phase portraits. Traveling waves corresponding to Lax 
shocks are given by trajectories for this system between a saddle point equilibrium 
and a second equilibrium, either a node or spiral. Oscillations are associated with a 
spiral equilibrium, at which the two eigenvalues are complex conjugates. Since either 
it + or u_ can be the saddle point, we examine eigenvalues for both, indicating whether 
there are oscillations behind or ahead of the wave. We find that the sign of p. is also 
significant. 

Intuitively, we expect that for small v > 0, there will be oscillations, since this 
is similar to underdamping of oscillations. For large v however, there should be no 
oscillations, corresponding to overdamping. We verify this distinction with a quanti¬ 
tative prediction of the boundary between the presence of oscillations in the traveling 
wave, and the absence of oscillations. The prediction is then observed in numerical 
simulations (see §6). 

If we linearize system (4.12) about an equilibrium u, we find eigenvalues A given 
by the formula 


A = 


1 

2 



2 

Consequently, the boundary between real and complex eigenvalues is ^- + 4 f u = 
0. Using the definition (4.12) of f(u,u-,s ) and the formula (4.3) for s, we obtain 
hyperbolas in the U -, u+ plane for each of u = U -, u = u+ : 

u = u_ : (2w_ + u + )(u_ — u + ) = — K- 

\ (4.13) 

u = u + : (2 u + + u_)(u + - u_) = -f^. 

These are represented in figure 4.3. The hyperbolas are oriented differently depending 
on the sign of /r. Also in the figure, we indicate the sectors corresponding to Lax shocks. 
There are several conclusions to be drawn from these figures. 

1. The dashed curves lie outside the sectors representing Lax shocks. Therefore, 
when n > 0, Lax shocks are non-oscillatory at U- and for /.t < 0, they are 
non-oscillatory at u+. 

2. For small enough v, Lax shocks are oscillatory at U- if fi > 0, and they are 
oscillatory at u+ if /.t < 0. 

Summarizing, we conclude that for v > 0 small enough, Lax shocks are oscillatory 
on the left (at w_) and non-oscillatory on the right (at u + ) if /z > 0, and the reverse is 
true if fi < 0. Therefore we can identify the shock orientation d = —sgn(/z), analogous 
to the orientation of a DSW (recall § 3.2). 

However, there is a significant difference between the location of the solid line 
hyperbolas in the two figures. In figure 4.3(a), the hyperbola crosses the line u+ = 
—u_/2, representing the shock in a shock-rarefaction wave, whereas in figure 4.3(b), 
the hyperbola is contained within the sector representing Lax shocks. Correspond¬ 
ingly, for fixed /.t > 0, and v > 0, rarefaction-shocks have an oscillation on the left, 
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whereas for /i < 0, rarefaction-shocks are monotonic. In both cases, weak shocks 
(those for which | u + — w_| is small) are monotonic. 

The superscripts in figure 4.2 are related to the appearance of oscillations for small 
6 > 0. For fi < 0, there are oscillations on the left of the wave, so that for it + < u_, the 
final oscillation has a maximum, indicating an S + wave, whereas in an S~ shock, we 
have u+ > it_, for which the final oscillation necessarily has a minimum. Note that 
in a solution designated SIS'" 1- , the shock is actually properly associated with the S + 
region with u + > w_. A similar remark applies to the SIS'” solutions. The S° shocks 
have monotonic traveling waves, in agreement with the remarks above concerning 
figure 4.3(b). 



Fig. 4.3. The hyperbolas delineating the boundary between oscillatory and non-oscillatory trav¬ 
eling waves. Lax shocks correspond to the shaded areas, (a) Region OSC where oscillations occur 
at U-, when ji < 0. (b) Region OSC where oscillations occur at «_)_ when // > 0. The dot on the 

u^-axis is at u_ = !v/2^/2|/If in (a), and at u_ = u/2y/Ji in (b). 

We have shown that there is an important distinction between Riemann problem 
solutions to the mKdVB equation depending on dispersion. When fi < 0, the effect of 
the dispersive term is the introduction of oscillations when S = v/-\/\fi\ is sufficiently 
small. Otherwise, the solutions are classical in the sense that shocks satisfy the 
Lax entropy condition. The classification diagram in figure 4.2(a) is independent 
of S > 0. While the ft > 0 case also yields oscillatory Lax shocks, this regime 
supports non-classical, undercompressive shocks that do not satisfy the Lax entropy 
condition. Here, dispersion plays a more significant role, which manifests as a in¬ 
dependent classification diagram (figure 4.2(b)). Note that the limit S —> oo (i.e., 
ft. —> 0, the purely diffusive case) of figure 4.2(b) leads to the same Riemann problem 
classification as in figure 4.2(a), i.e., the Riemann problem is independent of v when 

fi —¥ 0. 

5. Riemann problem for the mKdV equation. In this section, we present 
the classification of Riemann problem solutions for the mKdV equation, i.e., (1.4) with 
v = 0. Similar to the KdVB equation discussed in §3, the transition to the inviscid 
case v = 0 in the mKdVB equation is singular in the sense that although Lax shocks 
may have corresponding traveling waves for v > 0, they are not present for v = 0. In¬ 
stead, the Lax shocks are replaced by DSW solutions. In contrast to the KdVB/KdV 
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case, non-convexity of the hyperbolic flux in the case of the mKdVB/mKdV equations 
results in a much richer variety of solutions, which makes establishing the parallels 
between the diffusive-dispersive and purely dispersive cases rather non-trivial. One 
of the pertinent questions arising in this connection is: what is a conservative, purely 
dispersive counterpart of an undercompressive shock in the /i > 0 case? A sec¬ 
ond question concerns the identification of dispersive counterparts of composite wave 
structures such as shock-rarefactions, occurring as solutions of the mKdVB equation 
for both n > 0 and fi < 0 . 

The full classification of Riemann problem solutions for the Gardner equation, 
the extended version of the mKdV equation containing both quadratic and cubic 
nonlinearities, was constructed in [54]. In that paper, the Gardner equation was 
taken in the standard form w t ' + 6ww x i — 6 aw 2 w x > + w x ' x ' x ' = 0 , which can be 
reduced to the mKdV equation in the form (1.4) with v = 0, /jl = sgna by the change 
of variables u = y / 2|a|['u; — l/(2a)], t = t' , x = ~[x' — 3f'/(2a:)]sgna:. Here we 
adapt the classifications of [54] for both signs of a to the pure mKdV case (with the 
corresponding signs of fi) to enable detailed comparison with the counterpart mKdVB 
classifications outlined in §4. 

5.1. mKdV equation: DSW theory. To produce the classification of asymp¬ 
totic solutions to the Riemann problem for the mKdV equation, one needs first to con¬ 
struct elementary building blocks such as DSWs and rarefaction waves, and identify 
their admissibility conditions. For this, we first describe mKdV periodic solutions and 
their limiting configurations, and then develop modulation theory for mKdV DSWs. 
Unlike the well-established theory of diffusive-dispersive shocks for the mKdVB equa¬ 
tion, the counterpart mKdV theory is not directly available from a single source so 
we present here a unified account of previous results along with further development 
in the context of the theory of hyperbolic conservation laws. 


5.1.1. Traveling waves. Traveling wave solutions of the mKdV equation are 
available from a number of papers (see, e.g., [56, 75]) and can also be obtained as 
reductions of corresponding solutions to the Gardner equation [54]. 

Traveling wave solutions to the mKdV equation are sought in the form u = 
u(x — Ut), with speed U. Note that we have omitted the initial phase in the TW 
ansatz here as this does not contribute to the modulation solution of the Riemann 
problem (see Section 3.2 for the discussion of the DSW phase for the KdV case). 

After integrating twice, one arrives at the ODE (dropping the tildes) 

u 2 = sgn/z(u — ui)(u — u 2 ){u — u 3 )(u - U 4 ) = Q{u), (5.1) 


where 77 


x—Ut 


and the roots of the polynomial Q(u) satisfy the constraint 


4 

yi u.j = 0. 

i =1 


(5.2) 


While only three roots Uj are independent, it is still convenient to keep all four 
Uj's in the subsequent formulae to preserve symmetry of the expressions. In the 
modulationally stable case of our interest, all Uj’s are real. We assume the ordering 


Ui < U2 < M3 < W4 • 


The phase velocity U is expressed in terms of the Uj ’s as 

U = ~-(UiU 2 + UiU 3 + U1W4 + U 2 U 3 + U2U4 + M3W4). 


(5.3) 
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The ODE (5.1) can be viewed as a nonlinear undamped oscillator equation with the 
“potential” —Q(u). Since finite real-valued motion occurs only in the intervals of u, 
where Q > 0 (see figure 5.1), the structure of the solutions crucially depends on sgn/x. 




Fig. 5.1. Potential curve Q(u ) for traveling wave solutions of the mKdV equation: (a) fj, > 0; 
(b) fi < 0 

Integration of the ODE (5.1) for both orientations of the curve Q{u) shown in 
figure 5.1 leads generically to solutions expressed in terms of elliptic functions. The 
full classification of solutions arising when all four roots u 3 are real is presented in 
Appendix A. Some of these solutions are qualitatively analogous to the cnoidal wave 
solution (3.7) of the KdV equation exhibiting two distinct limits: a zero-amplitude 
harmonic wave and a soliton. At the same time, the fact that the potential Q{u) is a 
quartic polynomial gives rise to a number of qualitatively new features. For example, 
unlike the KdV traveling waves, periodic solutions of the mKdV equation exhibit 
solitons of both polarities for each sign of /. i . Moreover, there are entirely new families 
of solutions not encountered in the KdV case. For /i > 0 these are kinks - smooth 
heteroclinic transitions connecting two equilibria u\ = U2 and 1x3 = M4 (note that 
u\ A U 4 = 0 due to (5.2)), 

u = A — (r/,4 - Mi) tanh[— (m 4 - ui)rj\. (5.4) 

The traveling speed for both types of kinks is U = m| , which agrees with the classical 
shock speed (4.3). 

For /i < 0 there is a “non-KdV” family of nonlinear trigonometric solutions 
corresponding to the merger of two of the roots, say Mi and M2, while the solution 
oscillates between the two remaining troots M3 ^ M4, see (A.13). 

5.1.2. Modulation equations. The purpose of this section is to explore the 
structure of the mKdV modulation equations as a system of hyperbolic conservation 
laws. Modulation equations for the mKdV equation (both signs of \i) were first derived 
by Driscoll and O’Neil [18] using the original Whitham averaging procedure applied 
to the first three mKdV conservation laws. The resulting system describes the slow 
evolution of the roots Mi(x, t), u^ix, t),U3(x, t), uniquely defining the (locally) periodic 
solution ((A.l) for ^ > 0 or (A.11) for /z < 0) and its physical properties (period-mean, 
wavelength, amplitude etc.). By employing a number of non-trivial relationships 
between elliptic integrals it was shown in [18] that the linear combinations 

V = ^( u i + u^), A2 =-(«i + M3), A3 = 2 (m 2 A M3) 


(5.5) 
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are Riemann invariants of the mKdV modulation system. This result was inspired by 
the original finding of Whitham [97] for the KdV equation. It transpired later [82], 
[56], and more generally in [54] using the methods of finite-gap spectral theory [51], 
that the modulation system for the rnKdV equation can be mapped onto the KdV 
modulation system (3.9). Specifically, the mKdV modulation system has the form 


^+m(A)^ = 0, i = 1,2,3, (5.6) 

where the characteristic speeds HAj(A) < W 2 (A) < W\{\) are related to the KdV- 
Whitham characteristic speeds (3.10), (3.17) by 

M < 0 : Wi{ A) = V 4 -i(r), (5.7) 

M >0: Wi{\) = V 4 .i(r*), r* = (r 3 ,r 2 ,n). (5.8) 

Here, the KdV-Whitham Riemann invariants rj are expressed in terms of the mKdV- 
Whitham Riemann invariants \ :j as: 

M < 0: rq = 3A 3 , r 2 = 3A|, r 3 = 3A?, (5.9) 

fi > 0 : n = -3A?, r 2 = —3Al, r 3 = -3A§. (5.10) 


In each case ^ > 0,/x < 0, we have for the modulus mi,m 2 respectively (see (A.3), 
(A. 12)), the same formula 


= -, (5.11) 

r 3 - n 

which is identical to the KdV expression (3.8). The admissible sets of inverse formulae 
for Aj(r) (and hence for Uj( r)) are selected by the inequalities 

Ai < A 2 < A 3 , ri < r 2 < r 3 . 

Thus, for /x < 0 one has either 

Ai = -v^73, A 2 = -\AV3, A 3 = -Vnjs (5.12) 

or 

Ai = — \/r 3 /3, A 2 = —\J r 2 /3, A 3 = \J r 4 /3. (5.13) 

We note that relations (5.12) correspond to modulations of the periodic solution 
(A.11) occurring in the interval u 3 < u < 14 while relations (5.13) correspond to 
modulations of the counterpart solution occurring in iq < u < u 2 . 

Similarly, for /x > 0 one has either 

Ai = — ^i/3, A 2 = —\/—r 2 /3, A 3 = — \/—r3/3 (5-14) 


or 


Ai — — \/—t i/3, A 2 — — \/—T 2 / 3 , A 3 — r 3 /3. 


(5.15) 
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The asymmetry in the expressions (5.12) and (5.13) ((5.14) and (5.15)) is due to the 
use of three variables A j (5.5) replacing the four original quantities Uj (see (5.5)), 
related by the condition (5.2). 

Non-uniqueness in the determination of the vector A = (Ai,A2,As) for a given 
vector r = (r\,r2,T^) is due to the invariant transformation 

u —> — u, Ui~>—U 5 -i, i=l,2,3,4 (5.16) 


of the mKdV traveling wave solutions (see Appendix A) enabling two different sets 
of {uj} to be mapped onto the same set {r 7 }. Indeed, one can readily see that the 
sets (5.12) and (5.13) (as well as (5.14) and (5.15)) are related via the transformation 
(5.16). The choice of the actual set {Aj} is determined by the initial or boundary 
conditions for the Whitham equations. It is also worth mentioning that the mapping 
{Ai, A 2 , A 3 } 1 —t- {ri, r 2 , r 3 } can be viewed as the modulation theory counterpart of the 
Miura transformation connecting the KdV and mKdV equations [76]. 

The quadratic nature of the transformations (5.9), (5.10) implies that the mKdV- 
Whitham system (5.6) is neither genuinely nonlinear nor strictly hyperbolic. Indeed 
one can see from (5.7), (5.13) for /x < 0 and (5.8), (5.15) for /x > 0 that coinciding 
characteristic velocities are admissible: W^(Ai, —A 3 , A 3 ) = Ws(Ai, —A 3 , A 3 ), implying 
that system (5.6) is non-strictly hyperbolic. Non-genuine nonlinearity then follows 

[15]- 

It is instructive to derive the non-genuine nonlinearity property of the mKdV- 
Whitham system (5.6) directly, without invoking non-strict hyperbolicity. 

Let n < 0. Consider the derivative 


dWj 

d\j 



(5.17) 


Since ^ ^ 0, j = 1,2,3 due to genuine nonlinearity of the KdV-Whitham system, 
the derivative (5.17) can vanish only if the value A j = 0 is admissible for some j € 
{1,2,3}. It follows from (5.12), (5.13) that one has Xj = 0 only when ry = 0. From 
the inequality ri < r 2 < r 3 and relations (5.12), (5.13) one can see that xq = 0 
(i.e. A 3 = 0) is admissible. Therefore the mKdV-Whitham system is not genuinely 
nonlinear for /x < 0. The proof for /x > 0 is analogous and involves admissibility of 
r 3 = 0 in (5.14), (5.15). 

Corollary. The mKdV-Whitham system is genuinely nonlinear in any restricted 
A-domain not containing an open neighbourhood of A 3 = 0. 

One can see that non-convexity of the hyperbolic flux f(u ) = u 3 in the mKdV 
equation gives rise to complex structure of the associated modulation system, a feature 
not present in the classical, simple-wave DSW theory typified by the KdV DSWs 
described in §3.2. Indeed, the immediate consequence of non-convexity of f(u ) is 
non-genuine nonlinearity of the dispersionless limit equation (1.13), for which the 
characteristic velocity is V(u) = 3 u 2 and so the derivative V'(u) = 6 u vanishes 
at u = 0. This property is naturally inherited by the mKdV modulation system, 
simply because the dispersionless equation is an exact reduction of the modulation 
equations obtained both in the harmonic (to = 0 ) and soliton (in = 1 ) limits (see 
[22]). However, the mKdV modulation system (5.6) exhibits further complexity by 
possessing non-strict hyperbolicity. Non-strict hyperbolicity was also observed for 
the modulation systems associated with the fifth order KdV equation (the second 
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member of the integrable KdV hierarchy) [83], the complex modified rnKdV equation 
[61] and the Camassa-Holm equation [36]. For all these equations, as well as for the 
mKdV equation (1.4) with v = 0, n < 0, non-strict hyperbolicity is responsible for the 
occurrence of a new fundamental structure, a contact DSW , which will be described 
below. 

To conclude this section we note that, in spite of the same non-genuinely nonlinear 
inviscid/dispersionless limit shared by the rnKdVB and rnKdV equations, non-strict 
hyperbolicity does not enter the rnKdVB TW shock theory due to the scalar nature 
of the problem. 

5.1.3. Classical mKdV DSW solutions and admissibility conditions. 

The description of classical DSW solutions occurring in the Riemann problem for 
the mKdV equation requires obtaining self-similar simple-wave solutions of the mod¬ 
ulation system, in which all but one of the Riemann invariants are constant. Such 
solutions exist as long as the modulation system is genuinely nonlinear (see, e.g., 
[64]). The mKdV DSW admissibility conditions will then include, along with Lax- 
type causality conditions (3.16), an extra (convexity) condition ensuring that the 
modulation equations are genuinely nonlinear in the relevant restricted domain of de¬ 
pendent variables. Due to the hyperbolic, albeit not necessarily strictly hyperbolic, 
nature of the problem, the necessary and sufficient condition for genuine nonlinearity 
is that the interval defined by the initial discontinuity (1.12) must not contain an 
open neighborhood of the point u = 0. In that case, the mapping (5.9) between the 
KdV and mKdV modulations is one-to-one and the mKdV DSW modulation solution 
has the form of a self-similar expansion fan analogous to (3.15). To avoid ambiguity, 
we assume below that u + ^ 0 and consider the special case u + = 0 separately. Also, 
as we shall see, the case u_ = 0 does not involve DSW formation. Thus we shall be 
assuming that U-U+ ^ 0 in our formulation of the DSW admisibility conditions. 

For [i < 0 we have 

n=3u 2 + , r 3 = 3u 2 _, V 2 {3u 2 + ,r 2 ,3u 2 _) = |, (5.18) 

where the connection between r :j ’s and the mKdV modulation parameters A j is given 
by (5.12) or (5.13) (see the criteria below). The DSW edge speeds obtained from 
(5.18) are: 

S-=3u 2 + -3A, s+= 3 u 2 + + 2A, A = u 2 _-u\. (5.19) 

This modulation solution, inserted into the periodic TW (A.11) gives and approximate 
description of the mKdV DSW and is shown in figure 5.2(a) for (u-,u+) = (2,1). 
Similar to the KdV DSW case, the weak limit u(x, t, /i) as /x —> 0 does not coincide with 
the Lax shock (1.14) but is given by the period-mean value u(x, t ). Also, we re-iterate 
that the exact phase of the TW solution (A.11) in the mKdV DSW is undetermined 
within the leading-order modulation theory. This does not affect the macroscopic 
characterization of the DSW, which is fully determined via the parameters r-|, r 2 , r 3 
by the modulation solution (5.18). 

The solution (5.18) is subject to the admissibility conditions 

|u_| > |u+| , U-U+ > 0, (5.20) 

The meaning of inequalities (5.20) is elucidated when they are written in terms of the 
flux function f(u) = u , 

/'(«-) > /'(«+), />-)/>+) > 0. (5.21) 
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The first of the conditions (5.20) is a mKdV equivalent of the general DSW causality 
conditions (3.16), which guarantee that the dispersionless limit characteristics x = 
c±t + x o, where c± = f'(u±) = 3u±, transfer initial data into the DSW region. Note 
that this condition does not coincide with the counterpart Lax entropy condition (4.4) 
for the non-convex conservation law (4.1), which says that the wave is a Lax shock 
if and only if sgn(u_)«+ lies between — sgn(it_)u_/2 and sgn(u_)u_. This contrasts 
with the KdV case, where the DSW causality condition does coincide with the Lax 
entropy condition (see §3.2). 

The second condition in (5.20) and more generally in (5.21) is the convexity con¬ 
dition, which guarantees that the initial step range does not include the inflection 
point u = 0 of the hyperbolic flux, thus ensuring genuine nonlinearity of the modula¬ 
tion system for the solution involved. Therefore, the convexity condition guarantees 
a single-wave regularization of the initial step. One can see that for the KdV equa¬ 
tion, the second condition (5.21) is always satisfied so that any initial jump can be 
regularized by a classical DSW, provided the DSW causality condition holds. 

When sgn(u_) = sgn(w + ), this sign indicates which of the inverse mappings (5.12) 
or (5.13) should be used to obtain the modulation (Ai, A 2 , A 3 ) for the relevant periodic 
wave (solution (A.ll) for DSW - or its reflection via (A. 8 ) for DSW + ). If > 0 
then one needs to use relation (5.13) to describe modulations in the DSW + , while for 
U— < 0 one has to use (5.12) to describe modulations in the DSW - . The polarity p 
of the DSW ( see §3.2) satisfies 


P = ~sgn(/iu_), (5.22) 

where p = 1 corresponds to an elevation wave and p = — 1 corresponds to a depression 
wave at the DSW solitary wave edge (see §3.2). Formula (5.22) follows from the anal¬ 
ysis of the four admissible DSW configurations (two for each sign of /i) summarized 
in §5.2.1, 5.2.2 below. 

In more general terms, 

P = -sgn(/z/"(u_)). (5.23) 

The behaviors of characteristics for the modulation system in the DSW region are 
analogous to those for the KdV case as in figure 3.3. However, KdV exhibits only 
the DSW + for p < 0 (see figure 3.5). In contrast, (5.23) demonstrates that mKdV’s 
non-convex flux allows for both DSW + and DSW - regularizations without changing 
the dispersion sign. It is the product of the dispersion convexity and the hyperbolic 
flux convexity in (5.23) that determines a DSW’s polarity. On the other hand, the 
DSW orientation (see §3.2) is determined by the dispersion sign alone. For the mKdV 
equation with a fixed dispersion sign, DSW + and DSW - have the same orientation, 
i.e. the same relative positions (trailing, leading) of the solitary wave and harmonic 
edges. 

The mKdV DSW modulation solution for the case p > 0 is obtained from (5.18) 
by applying the transformation (3.17) and then using one of the inverse mappings 
(5.14), (5.15) to get the modulation (Ai,A 2 ,As) of the periodic solution (A.l). One 
uses the mapping (5.14) for a DSW + and (5.15) for a DSW - . The DSW admissibil¬ 
ity conditions have the same form (5.20) but the DSW polarity is opposite to that 
occurring in the case p < 0, see (5.22). The speeds of the trailing and leading edges 
for the DSWs of both polarities are (cf. (5.19)): s_ = 3+ A, s + = 3 u\ + 6A. 

In all cases, the leading (p < 0) or trailing (p > 0) soliton amplitude in the 
classical DSW is a = 2 (|w_| — |u + |). 
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Fig. 5.2. Dispersive shock waves for the mKdV equation, constructed at t = 1 from mKdV 
modulation theory, (a) Classical DSW, p = — 2-10 — 2 . (b) Non-Classical contact DSW, p = — 5-10 —4 . 
(c) Envelopes of the DSW (solid, martini glass shape, scaled for comparison) and the CDSW (dashed, 
bordeaux wine glass shape). 


If one or both of the admissibility conditions (5.20) fail for a given pair (u-,u+), 
then the regularization of a step (1.12) via a single classical DSW is not possible. If 
only the DSW causality condition is violated but the convexity condition holds true, 
then the dispersive resolution occurs via a rarefaction wave. The mKdV rarefaction 
wave solutions for |u + | > |u_| for t 1 have the leading order asymptotic form (4.11), 
same as for the mKdVB equation (1.4). In contrast to the counterpart rarefaction 
wave solution of the mKdVB equation, the corners of the mKdV rarefaction wave are 
smoothed by linear dispersive oscillations whose wavelength ~ \/\lA and amplitude 
decays like f -1 / 2 [66], similar to the asymptotic behavior in the KdV rarefaction wave 
described in Sec. 3.2. Note that only the n > 0 case is discussed in [66]. 

In the remaining two (/i < 0 and /i > 0) classical DSW-inadmissible cases involv¬ 
ing the breaking of the convexity condition, “non-classical” dispersive regularizations 
occur; these will be described in the following subsections. The borderline case when 
it + = 0 can be formally viewed as being of either classical or non-classical regulariza¬ 
tion type, and will be discussed briefly. 

5.1.4. Non-classical DSWs. We first consider the case p. > 0. Here, we already 
have an exact heteroclinic solution to the mKdV equation - the kink solution (5.4) 
- that violates the convexity condition. Indeed, the kink with the minus sign in 
(5.4) connects two states U- > 0 at x — > — oo and = — U- at a: -> oo so that 
U-U- 1 _ < 0. The antikink with the plus sign in (5.4) is the reflection of the kink under 
the transformation (A.8) and has the same non-convexity property. One can see that 
kinks violate both the DSW causality (5.20) and Lax entropy (4.2) conditions as well. 
Unlike the classical, expanding, DSW, the kink is a moving front approximating a 
shock wave, so must be evaluated against the Lax entropy condition. The Lax entropy 
condition agrees with the counterpart DSW causality condition (3.16) if one takes 

= s + = sk- However, the kink speed is sk = v?_ < f'(u± ) = 3u 2 = 3 u 2 , so that 
kinks are undercompressive as described in §(1.15). From the discussion of traveling 
waves corresponding to undercompressive shocks in § 5.1.1, we observe that the kink 
solution (A.9) is the strong <5 —0 limit of the diffusive-dispersive undercompressive 
shock solution (4.9), (4.10). Thus we identify kinks as undercompressive DSWs. Note 
that the undercompressive nature of mKdV kinks was discussed in [59]. In Sec. 5.2, 
we will incorporate them into the general context of the mKdV Riemann problem 
classification. 
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An exploration of solutions to the Riemann problem (1.12) for the mKdV equation 
with /i > 0 was undertaken in [11]. The amplitude of the DSW’s lead soliton was 
determined using an inverse scattering approach and the necessary and sufficient 
condition U-U+ < O-breaking of the convexity condition (5.20)-for the formation of a 
kink (termed a double-layer solution) was identified. Numerical simulations were used 
to identify kinks, classical DSWs, rarefaction waves, and double wave combination 
solutions. In the next section, we will provide the full classification of the rnKdV 
Riemann problem for both fi > 0 and /i < 0. 

The Riemann problem (1.12) for the rnKdV equation with /r > 0 and the particu¬ 
lar initial data u+ = — was studied in [ 66 ], where a detailed solution was obtained 
using matched asymptotic expansions. The long-time asymptotics of this solution is 
dominated by a kink, which agrees with the modulation theory approach where kinks 
are used as one of the ‘building blocks’ in the construction of long-time asymptotic 
solutions in more general Riemann problems, see [26], [54]. 

Now we turn to the case fj, < 0. Instead of kinks, a peculiar type of DSW occurs, 
which we term a contact DSW (CDSW). A CDSW, similar to a kink, connects two 
conjugate states > 0 at x —> —oo and u + = —U- at x —> oo (this is CDSW + ) 
but, unlike the kink, has an oscillatory structure resembling that of a classical DSW. 
Contact DSWs represent the dispersive counterparts of single-sided contact disconti¬ 
nuities (or sound shocks) known in the theory of hyperbolic conservation laws (see, 
e.g., [15]). 

The CDSW + solution for the mKdV equation is described by the special modu¬ 
lation A l 7 A 2 ,A 3 in which A 2 = A 3 so that the elliptic parameter m = 0 (see (5.11)) 
throughout the wave train, but the wave amplitude a = 114 — 113 = — 4A3 > 0. This 
sharply contrasts with KdV modulation theory where m = 
plitude. The CDSW + modulation solution has the form 

—Ai = U-, X2 = A3, 

W / 2(— U-, A3, A3) = Wz(—u~, A3, A3) = 

or explicitly, using (5.7), (3.12) and (5.9), 

Aa= ~;Wi +3 “ 2 -' 

Since A 2 = A 3 implies u\ = U 2 and the oscillations in the CDSW + are confined to 
the interval U 3 < u < U 4 , we conclude that the ‘carrier’ wave in a CDSW + is the 
nonlinear trigonometric solution (A. 13), and that is why this type of DSW has been 
termed a sinusoidal or trigonometric DSW elsewhere [75], [54]. At the leading edge 
of a CDSW + , Ai = A 2 = A 3 (i.e. u\ = U 2 = 113 ), which manifests a ‘bright’ algebraic 
soliton (A. 14) with amplitude a = —4Ai = 4 u- and speed U = 3Af. The algebraic 
soliton propagates on the background u = Ai which, according to the matching regu¬ 
larization conditions, must be matched to the constant state u = u + upstream. Then 
(5.24) yields that a CDSW + can only connect the conjugate states u = U- > 0 and 
u = u+ = —U- < 0. At the CDSW+ trailing edge, the wave amplitude a = 0, 
i.e., A 2 = A 3 = 0. Thus, the speeds of the CDSW + trailing and leading edges are: 
s_ = TV 2 ( — it— ,0,0) = — 2>u 2 _ , s+ = W 2 (—u~, — U-, — U-) = 3 u 2 _ = 3 u+. This mod¬ 
ulation solution, inserted into the periodic TW (A.13), is shown in figure 5.2(b). 
Similar to the classical KdV and mKdV DSWs, the exact phase of the CDSW is not 
determined by the modulation solution (5.24). 


0 implies zero wave am- 


x 

t 


(5.24) 



36 


G. A. EL, M. A. HOEFER, AND M. SHEARER 


The fundamental difference between the modulation solution (5.18) for a classical 
DSW and the solution (5.24) for a CDSW is that the solution of (5.24) represents 
a double characteristic fan, the corresponding distinct Riemann invariants being A 3 
and — A 2 . This solution is possible owing to the nonstrict hyperbolicity of the mKdV 
modulation system. Note that the trigonometric structure of the underlying mKdV 
traveling wave solution (A.13) is not the defining feature of CDSWs. See [61] where 
a qualitatively similar type of modulation occurs for a wave with m ^ 0. The more 
fundamental contact property is that the CDSW leading edge is a triple characteris¬ 
tic , tangential to the external dispersionless characteristic with characteristic velocity 
3m 2 . As a result, a CDSW + could be matched with a rarefaction wave at the leading 
edge. The composite CDSW-RW solutions will be considered in §5.2.2. We also men¬ 
tion the qualitative difference between the appearances of the envelopes for classical 
and contact DSWs: the classical DSW typically has a “martini-glass” shape due to 
the asymptotic behavior of the amplitude a ~ (x — X—) near the harmonic edge at 
x = X- = S-t, whereas the CDSW envelope defined by a = — 4 A 3 ~ y/x — X- has 
a “Bordeaux-glass” shape as shown in figure 5.2(c). We note that these envelope 
distinctions may not be obvious in direct numerical simulations due to contributions 
not included in the Wlritham multiple-scale solution, which could become important 
in the vicinity of the harmonic edge defining the appearance of the envelope. See [34] 
for the corresponding KdV analysis of higher order corrections to the harmonic edge. 

The description of a CDSW - , where the oscillations occur between the other pair 
of roots, Mi < u < M 2 , is analogous and involves the reflection of the traveling wave 
formulae (A.11) - (A.14) via the transformation (A. 8 ) and the relationships (5.13) for 
the modulation variables. 

5.2. Riemann problem classification. The Riemann problem classification 
for the mKdV equation on the m_-m+ plane of initial data is constructed by considering 
the admissibility conditions (5.20) and all the cases when one or both of them fail. As 
a result, the m_-m + plane is divided into eight equi-spaced sectors, each corresponding 
to a distinct fundamental wave regularization pattern (see figure 5.3). The separation 
lines are: m_ = 0, m+ = 0, u+ = m_ and u+ = — u_. Due to the symmetry inherent 
in the causality and convexity conditions (5.20), the boundaries between different 
solution types are symmetric under a 7 t /4 rotation. 

5.2.1. Classification for / 1 : > 0. The classification is shown in figure 5.3(a). 
We have already considered solutions corresponding to the classical DSW admissible 
regions defined by the conditions (5.20). The adjacent regions with |m_| < |m+|, where 
the DSW causality condition fails but the convexity condition holds, correspond to 
regularization via rarefaction waves described by (4.11). The line u+ = — u_ corre¬ 
sponds to a step regularization via a single kink (u_ > 0) or antikink (u_ < 0). These 
have also been described in §5.1.4. We now need to describe the remaining two cases 
corresponding to the sectors — u_ < u+ < 0 and m_|_ < — u_ < 0. Conjugate solutions 
corresponding to the opposite sectors are obtained via the reflection transformation 
(A. 8 ). 

Let — u_ < m+ < 0. The Riemann problem solution is a double-wave KjDSW + 
consisting of a kink and a DSW + connected by the constant state u = — u_. The 
kink speed sk = m 2 , and the speeds of the DSW + trailing and leading edges are 
s_ = m_ + 2u\ and s + = 6 u 2 __ — 3m^_ respectively. One can see that s_ > sk so the 
interval between the kink and DSW + increases with time. 

If U-j- < —m_ < 0, then the initial step regularization occurs via a kink and a 
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Fig. 5.3. Solutions of the Riemann problem for the mKdV equation. Legend: (R) - rarefaction 
wave; (DSW+) - ‘bright' DSW; (DSW~) - ‘dark' DSW; K - kink; (CDSW±) - ‘contact’ DSW ± ; 
(\) - intermediate constant state. 


rarefaction wave (formula (4.11) with the minus sign, and U- replaced with —u_) 
connected by a widening interval in which u = —u^. 

It remains to explain the singular transition between the patterns DSW - for 
0 < u + < U- and K|DSW + for — U- < u+ < 0, which occurs when crossing the 
horizontal axis u+ = 0. The trailing oscillation in the DSW - represents a depression 
(dark) soliton (A.6) having the amplitude 2(|u_| — \u+\) and propagating with the 
speed S- = 2u 2 + +u 2 _. When u + = 0 the soliton speed becomes equal to the kink speed 
Sk = u 2 _, and the left slope of the dark soliton becomes equivalent to a kink. The 
further decrease of u+ results in a ‘peeling’ off of the trailing kink from the remaining 
wave train, which becomes a DSW + . See [54] for a series of numerical simulations 
illustrating the analogous transition for the Gardner equation. 

5.2.2. Classification for ft < 0. The classification of the Riemann problem 
solution to the rnKdV equation with ft < 0 is shown in figure 5.3(b). The classical 
DSW regularizations corresponding to two opposite sectors of the U--U+ plane where 
the admissibility conditions (5.20) hold were described in §5.1.3. Similar to the ft > 0 
case, the adjacent regions with |u_| < |w + |, where the DSW causality condition fails 
but the convexity condition holds, correspond to regularization via rarefaction waves 
described by (4.11). The line u + = — u_ corresponds to a step regularization via a 
single CDSW+ (u_ > 0) or CDSW - (u_ < 0). These have also been described in 
§5.1.4. We now need to describe the remaining two cases corresponding to the sectors 
— U- < u + < 0 and it+ < — it_ < 0. The conjugate solutions corresponding to the 
opposite sectors are obtained via the reflection transformation (5.16). 

Let — U- < it+ < 0. The Riemann problem solution is a compound wave 
CDSW + DSW + consisting of a partial CDSW + and DSW + attached to each other. 
The line separating the CDSW + and the DSW + in the x-t plane is the character¬ 
istic s* = 6 u 2 + — 3 u 2 _. The trailing edge of the compound wave CDSW + DSW + is 
s_ = —3 v?_. The leading edge is s + = 2 u 2 _ + u 2 + . 

If u+ < —U- < 0, then initial step regularization occurs via a combination of a 
fully developed CDSW + connected at the leading edge to a rarefaction wave. 
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The classification of the regularization patterns in terms of the admissibility/non- 
admissibility conditions for both signs of /i is presented in Table 5.1. 


Causality (|u_| > « + |) 

Convexity (u-u+ > 0) 

Regularization /i > 0 {ft < 0) 

Yes 

Yes 

DSW" (DSW+) 

No 

Yes 

R (R) 

Yes 

No 

K|DSW+ (CDSW+|DSW+) 

No 

No 

K R (CDSW+ R) 


Table 5.1 

Classification of regularization patterns in terms of satisfaction/violation of the admissibility 
conditions (5.20) for U—U+ 7= 0. The DSW polarity is shown for u— > 0. 


6. Comparison. As can be seen by the very different analytical approaches 
outlined in §4 and §5, the Riemann problem solutions in the presence of dispersion 
with or without diffusion evolve in fundamentally different ways. However, the results 
do exhibit direct parallels that provide a useful bridge connecting diffusive-dispersive 
mKdVB dynamics to purely dispersive rnKdV dynamics. It is our aim in this section 
to identify the key differences and commonalities between these two models. 

6.1. Shock structure. The first and perhaps most obvious distinction between 
the two model types can be seen in the spatial structure of the solutions themselves. 
Figures 6.1 and 6.2 display numerical solutions to smoothed versions of the Riemann 
problem when /r = 1 and /.t = —1, respectively. See the Appendix for a description 
of the numerical method used. The first column of panels in each figure corresponds 
to purely dispersive dynamics with v — 0, whereas the second column corresponds to 
diffusive-dispersive dynamics in which v > 0. The text labels adjacent to each panel 
correspond to the classification in the (u_,u + ) regions of figures 4.2, 5.3. Due to the 
invariance of (1.4) under the transformation u —> —u, we restrict attention to the 
choice < u+, understanding that the analogous regimes when > u+ exhibit an 
amplitude reflection so + waves become _ waves and vice-versa. 

Many of the solutions exhibit oscillations due to the presence of dispersion. The 
presence of any diffusion damps the oscillation amplitudes. In contrast, the purely 
dispersive case leads to the continual generation of large amplitude oscillations show¬ 
casing the distinction between diffusive traveling waves approximating Lax shocks and 
expanding, dynamic DSWs. Even the rarefactions are subject to larger oscillation in 
the absence of diffusion. Notice, however, that both undercompressive shocks and 
kinks are monotone, exhibiting no oscillation. 

Due to the existence of two speeds, DSWs exhibit an orientation distinguishing the 
large amplitude, soliton edge from the small amplitude, harmonic wave edge. Another 
DSW characteristic is polarity, determined by the type of disturbance (elevation or 
depression) occurring at the soliton edge. Both orientation and polarity are inherited 
from the corresponding diffusive-dispersive wave when S = v/\f\fi\ is sufficiently small 
(recall §4.2), as shown in figures 6.1 and 6.2. The DSW orientation is determined by 
the sign of dispersion. As in the KdV equation, the only way to change the polarity of 
a DSW for a convex conservation law is to change the sign of dispersion, recall (5.23). 
A non-convex flux introduces the possibility of polarity change without changing the 
dispersion sign, as in figure 6.1(a). This can have physically important implications 
e.g. in internal wave dynamics (the distinction between dam-break and lock-exchange 
undular bores [26]) or in nonlinear optics/superfluids (the occurrence of bright solitons 
in a defocusing medium [14]). 
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t = 75 


t = 75 




a) fx = 1, v = 0 


b) fi = 1, v = 0.3 


Fig. 6.1. Numerical simulations of the Riemann problem for (a) mKdV and (b) mKdVB, with 
positive dispersion. The panel labels identify the types of waves, as in figures f.2 and 5.3. 


6.2. Characteristic diagrams. A useful means for understanding different 
wave solutions is through their characteristic diagrams. For scalar conservation laws, 
there is only one characteristic family. For the diffusive-dispersive case, this sin¬ 
gle characteristic family is sufficient to illustrate the TW solutions. However, in the 
purely dispersive case, shock solutions are described by a modulated, periodic TW via 
the Whitham equations with three distinct characteristic families. The characteristic 
description of single rarefactions, Lax shocks, and classical DSWs for the mKdV(B) 
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a) /x = — 1 , = 0 


b) fj, = —1, v = 0.3 


Fig. 6.2. Numerical simulations of the smoothed Riemann problem for (a) mKdV and (b) 
mKdVB, with negative dispersion. The panel labels identify the types of waves, as in figures f.2 and 


5.3. 


equation is similar to that for the KdV(B) equation shown in figure 3.3. 

Figures 6.3, 6.4, and 6.5 display characteristic diagrams and modulation solu¬ 
tions for double wave structures in the mKdV(B) Riemann problem classifications. 
The mKdV-Whitham characteristic families are denoted by Tj, i = 1,2,3 where is 
the set of characteristic curves x = x(t) satisfying = Wi, with W, the character¬ 
istic velocity in (5.6). Although we display generic diagrams, each characteristic is a 
legitimate one computed for a specific Riemann problem, i.e., these are not sketches. 
We now describe the implications of each of these diagrams. 

Figure 6.3 compares the shock-rarefaction, S + R in (a), and the contact DSW- 
rarefaction, denoted CDSW + R in (b-d), operable in the fi < 0 regime. Both the 
shock and the leading edge of the CDSW coincide with the trailing characteristic 
in the rarefaction, with speed s + . Moreover, the CDSW leading edge is a triple 
characteristic, i.e., s + = W\ = W 2 = W 3 . Here, all three modulation characteris¬ 
tics coincide, the reason for the terminology contact DSW. An additional feature of 
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Fig. 6.3. Characteristic diagrams for the composite shock-rarefaction (a) and contact DSW- 
rarefaction (b,c) solutions. The modulation solution (d) and associated three characteristic families 
Ti, r 2 , and T3 (b,c) are shown. The dashed curves correspond to T2 = T3- 


the CDSW distinguishing it from the classical DSW (recall figure 3.3) is the double 
characteristic family T 2 = T 3 , maintained across the entire contact dispersive shock 
structure. The remaining characteristic family passes straight through the CDSW, 
unchanged from the dispersionless characteristics with speed c_, so that c_ = s + also. 

Figure 6.4 depicts a CDSW + DSW + and allows for a side-by-side comparison 
of the non-classical CDSW and classical DSW characteristics. In this double wave 
structure, the CDSW modulation solution exhibits double characteristics W 2 = W 3 
up to its interface with the DSW where s* = W2 = W3. The Ti characteristic 
family carries information through the CDSW with speed c_ and into the DSW. 
The CDSW modulation solution (5.24) is unfolded by considering the alternative 
Riemann invariants A 3 < — A 2 < — Ai as in figures 6.3(d) and 6.4(d). This choice 
is natural for the CDSW + because its corresponding periodic TW solution (A. 13) 
satisfies U 3 < u < U 4 (recall figure 5.1). Our original choice of Riemann invariants 
(5.5) is more suited for TWs with u\ < u < U 2 - 

When > 0, composite waves involve undercompressive shocks, denoted E, for 
mKdVB, or kinks, denoted K, for mKdV. Their characteristics are shown in figure 6.5. 
In contrast to the fi < 0 double wave solutions, each wave in the composite solution 
for fi > 0 is separated by an intermediate constant state, the separation indicated 
by a vertical line |. The shock-rarefaction solution S°R for mKdVB, not displayed, is 
an exception in not having an intermediate constant; the speed of S° coincides with 
the characteristic speed c_ and forms the trailing edge of the rarefaction wave. A 
defining property of undercompressive shocks is the passage of a characteristic family 
through them. Figures 6.5(a,c) show how characteristics enter the undercompressive 
shock from the left and exit on the right with an upward deflection. In contrast, 
the corresponding kink solutions in figures 6.5(b,d) exhibit characteristics that pass 
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a ) b ) CDSW+DSW+ c ) 



Fig. 6.4. Characteristic diagrams for the three characteristic families Ti, r 2 , and T3 of the 
composite CDSW-DSW solution (a-c). The modulation solution is in (d). The dashed curves 
correspond to either r 2 = T 3 (between s~ and s*) or T 2 (between s* and s +). 


through with no deflection, just like one of the characteristic families for the CDSW 
solutions in figures 6.3(c) and 6.4(a). Through this similarity, the CDSW can also be 
thought of as undercompressive, however without a diffusive counterpart. Unlike the 
DSW and CDSW, the kink is a genuine, non-modulated TW solution of mKdV. Its 
effect in the modulation description of the K|DSW + of figure 6.5(e) is a discontinuous 
jump from U- to —u_, the intermediate constant connecting to the classical DSW. 
The KjDSW 1 * 1 is the only mKdV Riemann problem solution that incorporates both 
modulated (DSW) and non-modulated (kink) TWs. 

6.3. Zero Diffusion Limit. The limit v —> 0 + for solutions of the mKdVB 
equation is singular in that for every nonzero v , the long time shock behavior of the 
Riemann problem is resolved into a heteroclinic orbit or TW (4.9), (4.10). Diffusion 
introduces a relaxation mechanism to a steady configuration: TWs are attractors. 
But when v = 0, the only heteroclinic orbit available is the kink, valid only for very 
specific left and right states when /i > 0. So, the purely dispersive regularization for 
generic left and right states results in dynamic, oscillatory coherent structures. DSWs 
are unsteady, multiscale attractors that combine a periodic TW with an additional 
layer of self-similarity, the modulation solution. Nevertheless, we can see from the 
numerical results in figures 6.1 and 6.2 a correspondence between v > 0 and v = 0 
in most cases. It is therefore natural to introduce a mapping of the v > 0 mKdVB 
solutions to v = 0 mKdV solutions. As we will see, this mapping is multivalued and 
not one-to-one. 

We define the zero diffusion mapping for the Riemann problem by the classifica¬ 
tion figures 4.2 (v > 0) and 5.3 {v = 0) reproduced in figure 6.6. Common domains 
in the rt_-u+ plane demonstrate how diffusive-dispersive wave solutions map to the 
purely dispersive case. The mapping is summarized in table 6.1. The first and third 
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Fig. 6.5. Characteristic diagrams for Riemann problem solutions of mKdVB (a,c) and mKdV 
(b,d) with p > 0. The DSW modulation solution for Riemann invariants is shown in (e). 


negative dispersion fj, < 0 

positive dispersion ^l> 0 

R -)► R 

R —s> R 

S ± -► | 

DSW* 

or 

CDSW ± DSW ± 

[ DSW* 

S* -> l or 

{ K DSW* 

— 

SjS* -»• K|DSW± 

S*R-> 

r CDSW^R 

( or 

[ CDSW ± DSW ± 

£ R I f K R 

or >q; or 

S°R J [ K DSW* 


Table 6.1 

Zero diffusion limit mapping Riemann problem solutions of mKdVB to those of mKdV. 


quadrants of the U--U+ plane admit a single-valued, one-to-one mapping, comparable 
to the Riemann problem classifications for the KdV and KdVB equations with convex 
flux (recall figure 3.4). It is the remaining quadrants where non-convexity introduces 
novel features. The origin of the non coinciding boundaries for the Riemann problem 
classifications of mKdVB and mKdV can be traced back to the asymmetry intro¬ 
duced by the Rankine-Hugoniot and Lax entropy conditions (4.4) for mKdVB and 
the symmetric causality and convexity conditions (5.20) for mKdV. 

A direct mapping is possible when 0 < zt_ < u .|_ or U+ < U- < 0 for any fixed 
/i / 0. This implies that, as v —> 0 + , S* —> DSW*. Diffusive-dispersive Lax shocks 
map to DSWs in this regime. This provides an explanation for the similarity (within 
the regime considered) between the Lax entropy condition (4.2) and the DSW ad¬ 
missibility conditions (5.20), the new feature being the existence of two DSW speeds. 
The rarefactions are continuous and exhibit only weak discontinuities (i.e., in higher 
derivatives), hence their regularizations are the same to leading order. 
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Fig. 6.6. Replicas of figures f..2 and 5.3 with numerical solutions for interpretation of the zero 
diffusion mapping. 


In the case of positive dispersion /r > 0, the diffusive-dispersive classification 
approaches the same octant structure as the dispersive classification in the zero 
diffusion limit because <5 = —>• 0 + as v -» 0 + . This leads to the map¬ 

pings SIS* —^-KlDSW* and £|R —>K|R. Kinks are the purely dispersive analogue 
of diffusive-dispersive undercompressive shocks. This correspondence was similarly 
recognized in the context of two-layer fluids [59]. Numerical results in figure 6.1 show 
that the intermediate constant states associated with these double-wave structures 
attain different values according as v > 0 or v = 0. This is because a dispersive 
kink’s profile always admits the symmetry x —> —x, u —> — u. An undercompressive 
shock is not symmetric and the relation between the two constant states that it con¬ 
nects changes with the parameter S = v/ y yji, as described in §5.1.1. Nevertheless, the 
qualitative features between the diffusive-dispersive and dispersive solutions admit a 
natural correspondence. This correspondence with nonzero diffusion is neither one- 
to-one nor single-valued. There is the addition of the double-wave S°R, which does 
not exhibit an intermediate constant state nor does the shock exhibit any oscillation. 
Complicating the mapping further is the fact that the boundaries between the differ¬ 
ent wave structures when v > 0 and v = 0 do not line up. Then the mapping allows 
for E|R or S°R K|R or KIDSW*. 
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The negative dispersion fi < 0 case is complicated by the fact that, for any nonzero 
diffusion, the Riemann problem classification divides the U--u + plane into only six 
regions whereas the zero diffusion case involves eight. The new feature associated 
with zero diffusion is the existence of contact DSWs. As we have shown in previous 
sections, CDSWs have no “natural” diffusive-dispersive correlate. Nevertheless, the 
overlapping domains contained in the two octants 0 < it_ < — u+, —u+ < it_ < 0 of 
figure 6.6 imply that S ± R —»• CDSW ± R as v —» 0 + in this region. Shock-rarefactions 
become CDSW-rarefactions in the zero diffusion limit. A careful examination of 
the numerical results shows that the constant state connecting these double wave 
structures differ for the diffusive and non-diffusive cases. This is a manifestation of 
the difference between the Rankine-Hugoniot relation and the CDSW conditions. 

The zero diffusion mapping S ± —>CDSW ± DSW ± is implied by the regions 0 < 
u+ < —U -/2 and —«_/ 2 < u+ < 0 in figure 6.6. Lax shocks map either to contact 
DSWs when attached to a rarefaction or to a hybrid CDSW-DSW or to a pure DSW, 
depending on the relation between u- and u+. This implies that the zero diffusion 
mappings of S 1 * 1 and S ± R are multivalued. The top two panels in figure 6.1(a) show 
that, while there may be some subtle features introduced by CDSWs (recall figure 
5.2), generally it is difficult to distinguish them from pure DSWs. Recall that CDSWs 
arise due to nonstrict hyperbolicity at the level of the Whitham modulation equations, 
a feature of the purely dispersive case. Therefore it is perhaps not surprising that 
diffusive-dispersive shocks (TWs) do not discern such features. The introduction of 
a non-convex flux implies that Lax shocks map to either CDSWs or DSWs but their 
differences may be difficult to identify in practice. 

Aside from the ambiguity in the zero diffusion limit mapping of Lax shocks intro¬ 
duced by CDSWs when y < 0, the other limits described so far are natural. However, 
the remaining, unexplored regions of the U--U+ plane with regards to the zero diffu¬ 
sion mapping are u _/2 < —i t + < and U- < —w+ < u_/2 when /i < 0. Figure 6.6 
shows S ± R -A CDSW ± DSW ± in this region. But this mapping is peculiar because 
there is no rarefaction in the purely dispersive case. The resolution of this appar¬ 
ent contradiction can be understood by appealing to the small diffusion regime and 
relevant asymptotic time scales, which we undertake in the next subsection. 

6.4. Time scales and critical scalings. Motivated by the singular transition 
S ± R —A CDSW ± DSW ± just described, we now investigate the dynamics when the 
diffusion is small but nonzero. The existence of the parameter S = z^/VTa 4 ! i n the TW 
solutions of §4 provides a clue as to the critical scaling relationship between diffusion 
and dispersion. By “small diffusion”, we really mean 0 < v <C \f\lA- We now rescale 
space and time in (1.4) to demonstrate the time scales in which dispersion dominates 
the transient, intermediate asymptotic behavior, and the time scale on which the 
long-time behavior emerges, when diffusion and dispersion are in balance. In order to 
preserve the balance for transport on the left side of (1.4), we scale space and time by 
the same constant. Accordingly, let t = ar,x = ay. Then equation (1.4) transforms 
to 

u T + (u ) y = ~ u w + ^2 u yyv (6-1) 

In order for the dispersive term to balance the transport term, we set a = Then 

the diffusion coefficient of the first term on the right hand side becomes vI\/\y\- If we 
take this ratio to be small, then the diffusive term is negligible compared to dispersion 
over the time interval where t ~ \Z[/r[- 
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Fig. 6.7. Numerical simulation of the Riemann problem for U— = 1, = —0.75 comparing 

purely dispersive (left) and diffusive-dispersive (right) transient dynamics. 


Over long times, as the diffusive terms become significant, we expect solutions to 
converge to a combination of traveling waves and rarefactions, reflecting the combined 
effects of diffusion and dispersion. The time scale for such behavior emerges when the 
diffusive and dispersive terms are in balance. From (6.1) we see that the coefficients 
on the right hand side are equal when a = \ft\/v. The corresponding time scale is 
t ~ Since \ft\/v = ^ v\iAi we see that there is a transient regime 

■>J\ft\ -C t < \ft\/v over which the solution evolves from dispersion-dominated to the 
diffusive-dispersive behavior analyzed here. This transient evolution resembles the 
diffusion-dominated classical theory except in cases where undercompressive waves 
appear. 

As v approaches zero with fixed ft, the diffusive effects take longer to establish 
themselves, and in the limit v —> 0 + , the evolution is purely dispersive. The elongation 
of the transient dynamics time interval is suggestive of the transition between diffusive- 
dispersive solutions of the Riemann problem and purely dispersive solutions. Figure 
6.7 showcases this behavior for the numerical solution of a Riemann problem in the 
regime u_/2 < — u + < u- and ft < 0 where apparently S + R -A CDSW + DSW + as 
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v —► 0 + . We fixed fi = — 1 and evolved the same initial data, a smoothed transition 
from u_ = 1 to u + = —0.75 for v = 0 (figure 6.7(a)) and ^ = 0.01 (figure 6.7(b)). 
As predicted by the scaling analysis, the initial dynamics, e.g., t = 5 = 0{^/\p\), 
for both the dispersive and diffusive-dispersive models are essentially the same. The 
models initiate the formation of DSWs. But in the transient regime 1 = \f\ji\ <C1<C 
|/i|/i/ = 100, we observe the dynamics corresponding to a transition from a DSW to 
a shock-rarefaction. A large wave packet propagates backward, through the initiated, 
approximate DSW, leaving a steady diffusive-dispersive shock configuration adjacent 
to a rarefaction in its wake. Also predicted by the scaling analysis, at t = 100 = H/p, 
the diffusive-dispersive behavior has set in for the solution profile. This example 
indicates that the resolution of the singular limit v —> 0 + involves nontrivial, transitory 
dynamics. 

The transition S ± R —> CDSW ± DSW ± is a clear example of the singular behavior 
of the zero diffusion limit, demonstrating how to properly interpret the zero diffusion 
mapping in table 6.1 as a dynamic process. First, nonlinearity alone approximately 
describes the dynamics up to wavebreaking. Then, dispersion leads to the generation 
of oscillations and finally diffusion and dispersion together resolve the solution into 
a TW. As noted in §4.2, the resulting TW can be interpreted as an underdamped 
oscillator. The damping of this oscillator goes to zero with 6 , hence approaches a 
continually expanding dispersive shock wave. 

6.5. Small-diffusion regime and modulation equations. As we have seen, 
the analytical approaches used in the descriptions of diffusive-dispersive shocks and 
DSWs are quite different. The former is based on the analysis of an ODE for TW 
solutions of the mKdVB equation, while the latter employs nonlinear modulation 
theory, resulting in a system of homogeneous first order quasilinear PDEs. These 
two approaches can be reconciled in the case of small-diffusion 0 < v <C \/H- The 
reconciliation can be understood in the framework of an appropriate modification of 
Whitham modulation theory, in which the diffusive term in the mKdVB equation 
(1.4) is treated as a small perturbation of the mKdV equation. In this section, we 
briefly describe this approach for the well-studied case of the KdVB equation (3.4) 
and then outline its implications for the mKdVB case leaving a more careful analysis 
to future publications. 

In ‘perturbed’ modulation theory, the balance laws of the perturbed equation 
(e.g., the KdVB equation) are averaged over a family of periodic solutions for the 
unperturbed equation (i.e., the KdV traveling wave solutions (3.7) [30]; see also [20]). 
The result of this procedure is a non-homogeneous modulation system 

§ + q(r)|l= jS ,(r,, = 1,2,3, (6-2) 

where V 3 (r) are the KdV-Whitliam characteristic speeds (3.10) and the right-hand 
sides Rj{ r) arise from averaging the diffusive terms in the perturbed KdV balance 
laws. See [4], [41] for the explicit expressions in terms of complete elliptic integrals. 
The modulation system (6.2) for the KdVB equation with small diffusion can also be 
obtained using a direct multiple scales procedure, equivalent to Whitham averaging 
(see [77]). An effective general method for the derivation of modulation systems for 
perturbed integrable equations was developed in [52]. 

Unlike the unperturbed modulation system (3.9), the perturbed system (6.2) 
possesses TW solutions rj(x—st). It was shown in [4], [41] that the TW solution of the 
KdVB-Whitham system connects two constant, disparate states: u —» u± as x —> ±oo, 
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U- > u+, and satisfies the Rankine-Hugoniot condition s = ( u _ +u+)/ 2. Thus, this 
solution describes slow modulations of the classical diffusive-dispersive shock wave in 
the small diffusion regime. As a matter of fact, this outlined construction of the TW 
modulation solution is equivalent to a direct multiple scales perturbation analysis of 
the ODE (3.5) with S <C 1 for the TW solutions of the KdVB equation itself (see 
[49]). However, modulation theory provides a much broader platform, enabling one 
to describe different stages of diffusive-dispersive shock development (see [4, 20]). 

The general method [52] for the derivation of perturbed modulation equations 
via spectral finite-gap theory, and its applications to the rnKdV and Gardner equa¬ 
tions in [56], [54], suggest that the perturbed modulation systems for rnKdVB and 
KdVB are related by the same surjective mapping (5.9), (5.7) for fi < 0, and (5.10), 
(5.8) for ji > 0, as their unperturbed, purely dispersive, counterparts. This suggests 
that the modulation description of classical oscillating rnKdVB shocks can be derived 
from known solutions of the KdVB modulation system. However, the description of 
oscillatory contact shocks in shock-rarefaction complexes for /j, < 0 (see figure 6.7(b), 
t = 100) cannot be obtained from the KdVB modulation solutions and requires further 
analysis. 

7. Conclusions. The analysis of Riemann problems for the KdV(B) and 
mKdV(B) models reviewed here provides an explicit link between the theories of hy¬ 
perbolic conservation laws and nonlinear dispersive wavetrains. These are universal, 
asymptotic models of diffusive-dispersive Eulerian hydrodynamics. When the Eule- 
rian pressure law is non-convex, the modified KdV-Burgers equation emerges. The 
prototypical phenomena arising from the competition between diffusion, dispersion, 
and non-convex flux include non-classical wave features such as undercompressive 
shock waves, and double wave complexes such as shock-rarefactions. The purpose of 
this work has been to compare the singular, zero diffusion limits of these waves with 
their purely dispersive counterparts, identifying notable similarities and differences. 

An essential difference between diffusive-dispersive Lax shocks and classical dis¬ 
persive shock waves is that the former have fixed, traveling wave profiles whereas 
the latter waves have expanding, oscillating profiles. Consonant to this difference are 
the dissimilar shock jump conditions and their derivation. Nevertheless, comparing 
KdVB to KdV, there is a relatively simple, bijective zero diffusion mapping from Lax 
shocks and rarefactions with diffusion to DSWs and rarefactions without diffusion. 
The introduction of a non-convex flux in the rnKdVB and mKdV equations leads to 
a more complex, multivalued mapping. Now, the sign of the dispersion term plays 
an essential role. For negative dispersion and nonzero diffusion, the classical theory 
of conservation laws prevails. The new feature resulting from a non-convex flux is 
the shock-rarefaction, not essentially reliant upon the dispersive term in the equa¬ 
tion. In this regard, the small diffusion regime provides limited clues to the purely 
dispersive dynamics, where a new and unique feature arises, the contact DSW. These 
solutions occur due to the nonstrict hyperbolicity of the Wliitham modulation equa¬ 
tions. A CDSW is non-classical in the sense that either one characteristic family 
passes through it, as in a CDSW-DSW, or one characteristic family becomes tangen¬ 
tial to the characteristics of an adjacent wave, as in a CDSW-rarefaction. The former 
behavior is analogous to an undercompressive shock, which also has characteristics 
passing through, but occurs only in the positive dispersion regime. This behavior of 
characteristics in CDSWs contrasts sharply with the fact that all diffusive-dispersive 
shocks in the negative dispersion regime satisfy the Lax entropy conditions. As a re¬ 
sult, there is some ambiguity in the zero diffusion limit of Lax shocks. We explain and 
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resolve the ambiguity for small diffusion by identifying relevant time scales that in¬ 
clude a transient window during which the dispersive-dominated dynamics transition 
to a diffusive-dispersive balance. 

For positive dispersion, the diffusive-dispersive dynamics are fundamentally dif¬ 
ferent from the negative dispersion regime. Here, the role of dispersion is essential 
for the nonzero diffusion theory. This is due to the existence of non-classical TW 
solutions, namely undercompressive shocks, that are not available in the dispersion¬ 
less regime. The zero diffusion mapping identifies an undercompressive shock with a 
kink or undercompressive DSW. Kinks are monotone traveling waves satisfying the 
Rankine-Hugoniot shock conditions, and consequently do not exhibit the expansionary 
and oscillatory structure associated with classical DSWs. The Lax entropy conditions 
are not fulfilled for kinks because all characteristics pass directly through, without 
any deflection. 

This work emphasizes the universal scalar dynamics occurring in diffusive - dis¬ 
persive modifications (2.1) of non-convex Eulerian systems. Physically interesting 
dispersive and dispersive-diffusive systems exhibiting complex structure of solutions 
due to non-convex fluxes arise in fluid dynamics, plasma physics, magnetohydro¬ 
dynamics, nonlinear optics, Bose-Einstein condensates and magnetization dynamics. 
As examples of important bi-directional non-convex dispersive systems we mention 
Miyata-Choi-Camassa system for fully nonlinear internal waves [12], [26] and the 
cubic-quintic nonlinear Schrodinger equation arising in nonlinear optics and Bose- 
Einstein condensate dynamics [14]. 

The universality of the mathematical description developed in this paper is due to 
the generality of the rnKdVB equation (1.4), which is obtained via the well-established 
procedure of multiple scale expansions, and captures all essential properties of uni¬ 
directional, weakly nonlinear, long-wave diffusive-dispersive dynamics occurring in the 
neighborhood of the linear degeneracy point (2.9) for the background flow. While the 
rnKdVB description of the diffusive-dispersive dynamics of classical and non-classical 
shocks is based on the explicit calculations of TWs [47], the construction of the zero- 
diffusion DSW theory employs subtle integrability properties of the mKdV equation 
which are manifested, in particular, in the availability of Riemann invariants for the 
associated system of Whitham modulation equations. The extension of the present 
work to the fully nonlinear, vector (bi-directional propagation) case described by full 
non-convex Eulerian systems (2.1) or systems of similar structure will be a major 
development of the existing theory. 

We note that system (2.1) can be re-written using Lagrangian co-ordinates as a 
diffusive-dispersive modification of a general p-system (see, e.g., [15]). The theory of 
undercompressive shocks for non-convex p-systems was constructed in [68, 89]. The 
counterpart zero-diffusion theory is not available at present and its development will 
require an extension of the non-integrable methods of [22], [44] to the non-convex flux 
case. The inclusion of the possibility of non-convexity in the linear dispersion relation 
will introduce additional complexity to the Riemann problem classifications and their 
comparisons for diffusive-dispersive and pure dispersive cases. 

It is hoped that this work will inspire further examination and exploitation of 
the bridge between the well-established field of hyperbolic conservation laws and the 
growing field of dispersive hydrodynamics. 

Appendix A. Traveling wave solutions of the mKdV equation. 

We present the classification of the traveling wave solutions of the mKdV equation 
which are used in the construction of various DSWs in Section 5.1. These solutions are 
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obtained by integrating the ODE (5.1) for two cases /i > 0 and /x < 0. In describing 
traveling wave solutions we refer to two basic configurations of the potential curve 
Q(u ) shown in figure 5.1. 

If n > 0 (positive dispersion), then the oscillations occur in the interval «2 < u < 
m 3 , where Q(u) > 0 (see figure 5.1(a)), and the solution of (5.1) is expressed in terms 
of Jacobi elliptic functions as 


u = m 2 


(m 3 - M 2 )cn 2 (fl,TOi) 

1 — 113 - 112 sn 2 ( 6 >,mi) ’ 

U4—U2 v 7 - L/ 


(A.l) 


where 


9 = sj(us - mi)(m 4 - u 2 ) rj/2 , 


x — Ut 

VW\ 


(A.2) 


and 


m 1 


(m 3 - M 2 )(m 4 - Ml) 
(m 4 - M 2 )(m 3 - Ul)’ 


(A.3) 


is the modulus, 0 < m 1 < 1. The soliton limit m\ -A 1 can be achieved in one of two 
ways: when U 2 -A u\ or when u 3 —> 1 x 4 . 

When «2 ~A Mi, we obtain the “bright” soliton of elevation with amplitude a = 
us — Mi propagating against a constant background u = u \, 


m 3 - Ml 

M = Ml H- 4 - 4 — . 

cosh 2 0 - sinh 2 6 

U4-U1 

The soliton speed U is found from (5.3), (5.2) 

U = ^(3m 2 + 3MiM4 + m 2 ) . 

Analogously, for m 3 —a M 4 one obtains a “dark” soliton of depression 


(A.4) 


(A.5) 


U = U 4 — 


U4 — U2 

cosh 2 9 — sinh 2 6 ’ 

U4—U1 


(A. 6 ) 


having the amplitude a = M 4 — M 2 and propagating against the constant background 
u = U 4 with the speed 


U = ^(3m 2 + 3miM4 + m 2 ) . (A.7) 

Note that, due to invariance of the mKdV equation with respect to the reflection 
transformation u i-a —u formulas (A. 6 ), (A.7) can be obtained from (A.4), (A.5) by 
replacing 


M — A — M, Mi — A — Us—i , * = 1, 2, 3, 4. (A. 8 ) 

Thus the mKdV equation with fixed fj, > 0 admits the existence of solitons of both 
polarities. This is sharply distinct from the KdV dynamics where polarity of admissi¬ 
ble solitary waves is uniquely related to the dispersion sign (see §3.2). This property 
of the mKdV equation is due to the quartic nature of the potential curve Q{u) and 
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is ultimately related to non-convexity of the hydrodynamic flux of the dispersionless 
limit (1.13). 

If both U ‘2 —> ui and u 3 —> 114 then the polynomial Q(u) in the right-hand side 
of (5.1) has two double roots, which implies that the solution assumes the form of a 
kink/antikink, which can be obtained directly by integration of (5.1) with = 112 and 
u 3 = u 4 . Noting that due to (5.2) the double roots satisfy U\ + U 4 = 0 and choosing 
the constant of integration so that u = 0 at rj = 0 we obtain: 

u = ±-(«4 — u\) tanh [-(«4 — u\)rj\. (A.9) 

The lower sign corresponds to the kink with u —> 114 as rj —> —00 and u -A U\ at 
?? —> 00 ; the upper sign yields the “anti-kink” with u —> u\ as r/ —> —00 and u —> U 4 at 
?? —> 00 . The speed of kink propagation in both cases is found from (5.3) to be U = u\, 
which agrees with the classical shock speed (4.3), where u+ = ui, U- = U4 = —u\. 

We also present the small-amplitude asymptotics of the solution (A.l) when mi —> 
0 ( U 2 —> u 3 ). In this limit, the cnoidal wave (A.l) asymptotically transforms into a 
linear harmonic wave 

u = u 2 + \{u 3 ~ u 2 ) cos(kri), (u 3 - u 2 ) < u 2 , 

* _ 1 9 9 ( A - 10 ) 

k = \j(u 2 - Ui)(U 4 - U 2 ), u = ~{u\ + u\). 


Now we turn to the case of negative dispersion fi < 0. The potential curve 
configuration is shown in figure 5.1(b), and one can see that finite, real-valued motion 
can occur either in the interval u\ < u < 112 or m u 3 < u < U4. We only consider 
here the second case; the results for the first one can be obtained by applying the 
transformation (A.8). 

Let u 3 < u < U4. Integration of (5.1) yields 


(u 4 - u 3 )cn 2 (6»,m 2 ) 

U — Wq -, 

l+i^fsn 2 (0,m 2 )’ 

where 6 is defined by (A.2) and the modulus (cf. (A.3)) 

_ (U4 - U 3 )(U2 - Ml) 

(■U4 - u 2 )(u 3 - Ui)' 

In the soliton limit U 3 112 (m 2 —> 1), we get 


U4 - u 2 

u = U2 H-4- k— . 

cosh 2 0 + sinh 2 9 

U2—U1 


(A.ll) 


(A.12) 


This is a “bright”, elevation soliton with amplitude a = 114 — U2- Its propagation 
speed is found from (5.3), (5.2) to be 


U — — ( u 2 + 2U2U4 + Su 2 ). 


The limit m 2 —> 0 can be reached in two ways. 
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(1) If it 3 —> it 4 , we get asymptotically 


u = u 3 + -{u 4 - u 3 ) cos(kri), it 4 - u 3 < u 3 , 
k = \J(u 3 - ui)(u 3 - tt 2 ), U = ^(3u| + 2 m 2 U 4 + u\). 

This is a small-amplitude harmonic limit analogous to (A.10) obtained for [i > 0. 

(2) If U 2 = ui, but it 3 7 ^ U 4 then we arrive at the nonlinear trigonometric solution 


where 


u = u 3 


U4 - u 3 


1 + 


U4—U1 

U3-W1 


tan 2 6 ' 


8 = \J (zt 3 - iti)(u 4 - tti) 77 / 2 , U = -(3«3 + 3u 2 + 2 U 3 U 4 ), 

O 


(A.13) 


and the amplitude a = U 4 — u 3 is generally order unity. The nonlinear trigonometric 
solution (A.13) has no analogs in the KdV theory. Indeed, the very existence of the 
nonlinear trigonometric solution is possible due to the quartic nature of the potential 
function Q(u), i.e. due to the non-convex hyperbolic flux of the rnKdV equation. 

If u 3 —> U2 — U\, then 6 —» 0 and solution (A.13) transforms into an algebraic 
soliton of elevation, 


u = u 1 


U4 U\ j. _ 2 

1 + (u 4 -ui)V/4’ ~~ Ul ' 


(A.14) 


Note that the algebraic soliton propagates with the characteristic speed of the hy¬ 
perbolic limit (1.13) evaluated at the background u = u 3 . The counterpart “dark” 
solutions corresponding to the interval iq < u < u-i are obtained from (A.ll) - (A.14) 
by applying the reflection transformation (A.8). 

Finally, we mention a family of periodic solutions corresponding to the potential 
curve configuration where the roots u 3 and u 3 are complex conjugates. Such solutions 
are modulationally unstable [18] but they do not arise in the Riemann problem so we 
do not describe these solutions here. See [25] for the proof of hyperbolicity of the zero 
dispersion limit of the rnKdV equation with /a < 0. 


Appendix B. Numerical method. 

Here we describe the numerical method of solution for the KdV(B) and rnKdV(B) 
equations. The stiffness of the third order derivative term, the long time integration 
required, and the large domain needed to properly resolve the wave solutions make 
this a challenging computational problem. A high order, integrating factor (IF), 
pseudospectral Fourier method that incorporates differing, constant boundary values 
is developed, motivated by Trefethen’s approach to solving KdV [92]. The IF-RK4 
method, utilizing the fourth-order Runge-Kutta algorithm (RK4), has been shown 
to perform only slightly slower than other state-of-the-art methods for the KdV and 
Burgers equations [57]. Its ease of formulation and implementation make IF-RK4 
attractive. 

Interest here is in the long time behavior of the Riemann Problem 


Ut + c 0 u p u x = vu xx + \JUU XXX , i£l, t > 0, 


(B.l) 


u(x, 0) 


u- x < 0, 
u+ x > 0. 


(B.2) 
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The parameters include Co £ {1,3}, p £ {1,2}, v > 0, p £ R. 

We consider a truncated domain x £ [—L, L] and assume that, to an appropriate 
level of approximation, u(±L,t) = u± and d x u(x,t) is smooth and compactly sup¬ 
ported on [—L,L\ for each t £ [0,T], some T > 0. Then, the function v(x,t) defined 
by 

v(x,t) = u x (x,t), u(x,t) = J v(y,t)dy + u-, (B.3) 

is smooth and compactly supported in this domain. Differentiating (B.l) gives a 
nonlocal equation for v 


v t + c 0 (u p v) x = vv xx + nv xxx . 


(B.4) 


Since v(x,t) is localized within [— L, L], we can treat its periodic extension with the 
Fourier series expansion 


v(x,t) = ^2v n {t)e lknX , v n (t) 

n 


l 

2L 


j(x,t)e 


—ik n 


mr/L. 


The Fourier coefficients v n satisfy the differential equations 


Letting 


dv n 

d t 


+ ipk^Vn 


-ic 0 k n (v?v) - vk 2 n v n 


V n (t) = e 


— eivKti; 


i(t), 


then 


dV n 
d t 


-ic 0 k n e ipk ™ t tuPv) n - uk 2 V„. 


(B.5) 


A suitable truncation of the Fourier series is chosen, \n\ < N, and the corresponding 
eqs. (B.5) are then stepped forward in time by the standard RK4 method. The 
stiffest, third derivative term from (B.4) is absorbed into the transformation to V n , 
thus allowing for stable timestepping with larger At. We do not incorporate the 
second derivative term into the phase transformation because it leads to exponentially 
large/small terms during the calculation and causes numerical instability. 

A calculation shows that the accumulation function in (B.3) satisfies 


u{x,t) = ^2 

n^O 


Vn{t) c ik n x 

ik n 


1 

2L 



xv(x, t) dx + {x + L) 


u+ — u 
2L 


+ U-. (B.6) 


Equation (B.6) is suitable for numerical computation and is used for the reconstruction 
of u(x, t) where all Fourier series coefficients are computed using the FFT. The integral 
in (B.6) is computed with the trapezoidal rule. So long as the solution v(x,t) remains 
localized within [—L, L], this method is spectrally accurate. 

The initial condition (B.2) is smoothed by a hyperbolic tangent profile with width 

one. 

For the computations presented in this work, the numerical method was thor¬ 
oughly tested and validated using known, exact traveling wave solutions. Generally, 
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the purely dispersive computations required more care due to the radiation of small 
amplitude dispersive waves to the boundary. We monitored the conserved quantity 
f_ L v(x,t) dx to be within 10 -13 its nominal value u+ — U- across all simulations. 
The deviation of u at the boundaries from u± was maintained below 10 -4 (usually 
well below that). The numerical parameters were informed by the validation studies 
and chosen based upon the particular Riemann problem considered. Typically, we 
used Ax e {0.05,0.1,0.2}, At e {1CT 3 ,1CT 4 ,5 • 1CT 5 }, and L € {400,800, 1600}. 
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